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Abstract. We analyse the aging dynamics of the one-dimensional Fredrickson- 
Andersen (FA) model in the noncquilibrium regime following a low temperature 
quench. Relaxation then effectively proceeds via diffusion limited pair coagulation 
(DLPC) of mobility excitations. By employing a familiar stochastic similarity 
transformation, we map exact results from the free fermion case of diffusion limited 
pair annihilation to DLPC. Crucially, we are able to adapt the mapping technique 
to averages involving multiple time quantities. This relies on knowledge of the 
explicit form of the evolution operators involved. Exact results are obtained for 
two-time correlation and response functions in the free fermion DLPC process. The 
corresponding long-time scaling forms apply to a wider class of DLPC processes, 
including the FA model. We are thus able to exactly characterise the violations of the 
fluctuation-dissipation theorem (FDT) in the aging regime of the FA model. We find 
nontrivial scaling forms for the fluctuation-dissipation ratio (FDR) X = X(t w /t), but 
with a negative asymptotic value X°° = —3ir/ (6tt— 16) w —3.307. While this prevents 
a thermodynamic interpretation in terms of an effective temperature, it is a direct 
consequence of probing FDT with observables that couple to activated dynamics. The 
existence of negative FDRs should therefore be a widespread feature in non mean-field 
systems. 
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R(t,t w )=T 



(2) 

h=0 



Introduction 

A generic feature of glass-formers is the rapid increase of their relaxation time with 
decreasing temperature. When quenched below the laboratory glass transition, where 
relaxation times exceed experimental timescales, these systems fall out of equilibrium. 
They change from equilibrated fluids to nonequilibrium amorphous solids [lj. The 
"waiting time" t w elapsed since preparation of the state then sets a timescale for the 
subsequent relaxation: the system ages [2] . A full understanding of this nonequilibrium 
phenomenon remains a central theoretical challenge. 

Mean-field models of structural and spin glasses [31 H], in which relaxation times 
strictly diverge at the glass transition, have delivered important insights into the 
dynamics of aging. Most notably, these systems satisfy a generalized fluctuation- 
dissipation theorem (FDT) in their nonequilibrium, aging state. This is defined in 
terms of the two-time connected correlation function for a generic observable A, 

C(t,t w ) = (A(t)A(t w )) - (A(t))(A(t w )), (1) 

with t >t w , and the corresponding two-time response function 

S(A(t)) 
Sh{t w ) 

Here h denotes the thermodynamically conjugate field to the observable A so that 
the perturbation of the Hamiltonian is 57i = —hA. Note that we have absorbed the 
temperature T in the definition of the response. The associated generalized FDT is then 

R(t, t w ) = X(t, t w )—C(t, t w ), (3) 

with X(t,t w ) the fluctuation-dissipation ratio (FDR). At equilibrium, correlation and 
response functions are time translation invariant, depending only on At = t — t w , and 
equilibrium FDT imposes that X — 1. This is no longer true in nonequilibrium systems. 
But the definition of an FDR through Eq. ([3]) becomes nontrivial for aging systems: 
in the mean-field spin glass models 0, H] its dependence on both time arguments is 
only through the correlation function X(t,t w ) ~ X(C(t,t w )) at large times. This led 
to the introduction of timescale dependent effective temperatures and a possible 
thermodynamic interpretation of aging j3jlU[6]. 

In many systems of physical interest, however, such as liquids quenched below 
the glass transition or domain growth in disordered magnets [7], the dynamics is 
not of mean-field type. Crucial new features are activated processes and spatial 
heterogeneity 0, [9] [10] [IT] . Some experiments and simulations [12] have nonetheless 
seemed to detect a mean-field aging regime. On the other hand, theoretical studies 
have found ill-defined FDRs [13], [H] , non-monotonic response functions [T5], [16], (TT] [18] , 
observable dependence [T91I20]. nontrivial FDRs without thermodynamic transitions [21] 
[22] and a subtle interplay between growing dynamical correlation lengthscales and FDT 
violations [231 [21]. Also in experiments deviations from mean- field expectations have 
been observed, with for example anomalously large FDT violations associated with 
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intermittent dynamics [25J. It thus remains an important task to delineate when the 
mean-field concept of an FDR-related effective temperature remains viable. 

In this paper, we focus on the one- dimensional Fredrickson-Andersen (FA) model 
[26| [27]. Its relaxation time follows an Arrhenius law at low temperatures, and it is 
the simplest kinetically constrained model for glassy systems which displays dynamical 
heterogeneity [28l [29], [30l [3T] . We use it to study systematically the impact of activated, 
and in this sense strongly non mean-field, dynamics on FDRs and associated effective 
temperatures. In addition, the FA model dynamics become critical at low temperatures 
[32| 1331 [3l], with a diverging dynamical lengthscale. Our work is also of relevance, 
therefore, to the study of FDRs in nonequilibrium critical dynamics. Here, FDRs 
play the role of universal amplitude ratios, which makes them important markers for 
distinguishing dynamic universality classes [35] . 

A first numerical study of FDT violation in the one-dimensional FA model appeared 
in [13] and suggested ill defined FDRs. The authors used disconnected correlations, 
however, which compromises their results since FDT would not be recovered in 
equilibrium. More recent numerical studies [211 ES] based on the correct, connected 
correlations indicate a very different picture: no detectable violation of the equilibrium 
FDT is observed in spite of the nonequilibrium aging dynamics. We show below that this 
is a consequence of the specific scaling in the quasi-equilibrium regime [37J, which acts 
to obscure genuine aging contributions. These contributions are nevertheless present 
as we will show, leading to FDT violations with well-defined FDRs in the activated 
aging regime of the FA model. We take advantage of observable dependence of the 
FDR in order to gain more direct access to FDT violations, by considering global rather 
than local observables. Our discussion also elucidates the physical origin of negative 
dynamical response functions, and predicts the generic existence of negative FDRs for 
observables directly coupled to activated processes. A partial summary of our results 
can be found in [38] . 

Our analysis of aging in the FA model is based on a close connection with diffusion 
limited pair coagulation (DLPC) processes. The research field of diffusion limited 
reaction processes is well developed [39] and offers powerful analytical tools [101 W\\ . 
These include, in particular, stochastic similarity transformations [42l [43] . We use the 
latter to map exact results for the free fermion case of diffusion limited pair annihilation 
(DLPA) to DLPC. In particular, we show how the mapping technique can be adapted 
to averages involving multiple time quantities. This relies on knowledge of the explicit 
form of the evolution operators involved. We derive simple and explicit propagators for 
certain observables in DLPA and DLPC. In the free fermion DLPC process this yields 
new and exact results, for instance, the two-time particle correlation functions. Our 
propagators are of generic value for the analysis of DLPA and DLPC, however, and 
could serve in further studies of these processes. 

The paper is organized as follows: in Sec. [I] we recall the definition of the FA model 
and introduce the matrix formalism used throughout this manuscript. A discussion 
of the dynamics in the different stages of relaxation is given, and we establish the 
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connection to DLPC. The mapping between DLPC and DLPA is then stated in Sec. [2J 
We introduce suitable sets of observables, which are related in a simple manner by 
the mapping, and give expressions for their dynamics in the free fermion processes. 
Arguments are presented that establish the broader validity of the predictions for the 
long-time scaling behaviour. The generalisation of the mapping procedure to multiple 
time averages and propagators is presented in Sec. [3j Using the DLPC propagator, 
expressions for two-time correlation and response functions are derived in Sec. 13.11 and 
Sec. 13.21 Response functions in the FA model comprise two contributions representing 
different physical effects and correspondingly Sec. 13.21 is split into two subsections. 
Readers more interested in results than in their derivation may skip Sec. [2] and Sec. [3] on 
first reading, except for the beginning of Sec. 13.21 where our decomposition of response 
functions is introduced. Results for FDT violation in the activated aging regime of 
the FA model are then presented in Sec. H] and Sec. [5] and deal with local and global 
observables, respectively. Each of these sections comprises three subsections, discussing 
separately the dynamics of correlations, response functions and the resulting FDRs. We 
conclude with a summary and discussion of our results in Sec. El Two short appendices 
contain useful mathematical and technical material. 

1. Matrix Formalism 

In this section we summarize the background for our subsequent analysis of the aging 
dynamics in the FA model. After recalling the definition of the FA model we review the 
standard operator formalism widely used in the literature for studying reaction diffusion 
problems. We then translate the FA model into this representation and discuss its 
effective dynamics in the different stages of relaxation after a quench. 

The one-dimensional FA model [26J is defined on a linear lattice of size N subject 
to periodic boundary conditions. Assigned to each site i is a binary variable n.j e {0, 1}, 
with rii = 1 representing a mobile fluid region and = an immobile one. The state 
of the system is characterised by n = (ni, n 2 , . . . , n N ). The probability p(n, t) of being 
in state n at time t evolves according to the generic master equation 



where d t denotes the time derivative and w(n — > n') the rate for transitions from state 
n to n'. In the FA model these transition rates are given by 



The Kronecker deltas, S nin > = 1 if n = n' and zero otherwise, express that transitions 
occur only between states n,n' differing by a single flip — > 1 — rij. Under Glauber 
dynamics at temperature T the (unconstrained) flip-rates are 






(5) 




(6) 
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and obey detailed balance with respect to the noninteracting Hamiltonian 7i = Yli n i- 
Here e = (n^) = 1/(1 + e 1//T ) is the equilibrium density of mobility excitations; at low 
temperatures this is small and so the excitations can also be thought of as defects. The 
key ingredient of the FA model is the kinetic constraint 

fi(n) = Ui-i + n i+1 . (7) 

The dynamics on each site i require facilitation by an adjacent mobility excitation 
7ii_i = 1 or n i+ i = 1. Since /»(ra) is independent of it preserves detailed balance; 
the equilibrium state of the FA model is trivial. However, the kinetic constraint 
induces nontrivial relaxation dynamics with a rich phenomenology including dynamical 
heterogeneity [28l EHJ [30l EH [32l [33l EH SH SHI [46] . 

The discussion throughout this manuscript is based on a standard operator forma- 
lism. In the following we briefly summarize the main ideas; more details may be found 
in [4TJI |4"3] . One associates with each site % two orthonormal states |0) and |1), or \rii) 
for short, representing the mobility rij. Then an orthonormal basis for the configuration 
space of the entire system is 

\ n ) = K) with (n\n') = Y[S niiTl! .. (8) 

i i 

Probabilistic states \P(t)) are defined by 

|P(f)) = 5>(tM)|n), (9) 

n 

from which the probability for any particular configuration n is extracted via p(n, t) = 
(n\P(t)). One now multiplies the master equation Eq. by \n) and takes a sum over 
n. The result may be cast in the form 

d t \P(t)) = W\P(t)), (10) 

with W the master operator 

W = ^2w(n -> ri) (\n')(n\ - \n)(n\). (11) 

n,n' 

Conservation of probability is compactly expressed by (e\P(t)) = 1 if one introduces 

(e| = ^H = ((0| + (l|f Jv . (12) 

n 

This is the bra ground state of W because (e\W = from Eq. (flOl) . The corresponding 
ket ground state is by definition the equilibrium state |P eq ) so that W is generally a 
non-hermitian operator. The formal solution of Eq. ( flOl) for initial state |P(0)) is 

\P(t)) = e wt \P(0)), (13) 

with e wt the evolution operator. To calculate the expectation value of an observable 
A(n) one introduces the associated operator A = X]n yZ K n )l n )( n l in terms of which 

(A(t)} = (e\A\P(t)}. (14) 
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It is straightforward to verify that indeed (e\A\P(t)) = ^ n A(n)p(n,t). This completes 
our recap of the operator formalism. 

The FA master operator W FA is obtained by substituting the transition rates Eq. ([5]) 
into Eq. (TIT]) . To state it explicitly it is useful to introduce the number operator 

J2n t \n)(n\ = (\l)(l\) t . (15) 



ri; 



The notation (X)i indicates that the operator X lives in the subspace of site i, that is, 
(X)i = l®* 1 ' -1 ) ® X ® with 1 = 1 0> <0 1 + |1)(1| the diagonal operator; we will not 

distinguish the identity operator 1 from the number 1 in our notation as the meaning 
is always clear from the context. We also introduce the flip operator 

fl = (|l><0| + |0)<l|) 4 , (16) 

which maps states \n) onto F^n) = \rii, . . . , 1 — rii, n i+ i, . . . , n^). It is important 
to note that Ft and hi do not commute. But of course [Fi, hj] = if i ^ j since both 
operators act locally. Other obvious and useful properties are hf = hi and Ff = 1. The 
FA master operator now reads 

w fa = - m^, (17) 

i 

where fi = h^i +h i+ i and Wi = e(l — hi) + (1 — e)hi. As will become clear immediately 
it is useful to substitute the explicit form of f and, through a shift of the summation 
variable, rewrite the master operator as 

W FA = A > W f K = - !) W+i + {Fi - l)h i+ iWi. (18) 

i 

The key feature of Eq. ffTSl is that Wf A operates only in the subspace \ni,n i+ i). It 
therefore represents a pair reaction process. In the analysis of pair reaction processes it 
is conventional to choose an explicit representation for the basis |0) = Q and |1) = (°). 
All operators then have appealingly simple matrix representations. From Eqs. ( fl5l) and 
(|T6|) one has hi = (J J), and F { = (° ).- In particular, the FA master operator in the 
\ni,n i+ i) subspace becomes 

/ \ 

0-/3 7 



FA 



(19) 



-(3 7 

with (3 = e and 7 = 1 — e. The rows and columns of this matrix correspond to the 
configurations 00, 01, 10 and 11 from top to bottom and left to right, respectively. 
Consequently (3 is the branching rate, i.e. the rate for processes 01 — ^ 11 and 10 — ^ 11, 
while 7 is the rate for coagulations 11 — > 01 and 11 — > 10. 

We emphasize that the transitions 101 — > 111 or 111 — > 101 are accounted for 
correctly within this pair reaction description. In the case 101 — > 111, for instance, 
one has f(n) = 2 and Wi(n) = e according to Eqs. and ([7]), which amounts to a 
transition rate of 2e. Equivalently, from Eqs. (fl8l) and ( 1191) the process 101 — > 111 is 
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achieved by either branching to the right 101 — > 11 1 or to the left 101 — > 1 11 so that its 
total rate is again £ + £ = 2£ as required. The pair reaction picture is of course exact. 

We note that a modified version of the FA model with kinetic constraint fi(n) = 
7ij_i + n i+ i — ni_in i+ i has often been used in the literature since it was first proposed 
[17] . It has certain advantages for simulations because f G {0,1} is just a binary 
variable indicating whether the constraint of having at least one mobile neighbour is 
satisfied. However, for analytical work the presence of the quadratic term ni_in i+ i is 
awkward and prevents in particular a pair reaction representation like the one above. 
Nevertheless the scaling behaviour in the modified FA model is expected to match with 
the original version [26J discussed here. 



1.1. Effective Dynamics 

The branching and coagulation process Eq. ( TT9I) gives rise to the remarkably rich 
dynamics of the FA model. Here we are interested in the evolution after a quench from 
a random uncorrelated initial state (equilibrium at T = oo) to some low temperature 
T <^ 1 where e ~ e -1 / T < 1. It is well known that the dynamics then evolves in distinct 
stages with separated time scales [27] . 

On the 0(1) time scale activated states relax into blocked configurations. That is, 
the dynamics is controlled by coagulation processes 11 — > 01 or 10 which proceed with 
rate 7 ~ 1 until no further activated states 11 remain. This process has been analyzed in 
[48. 49]. The density of excitations drops from its initial value of c — | to c exp(— cq). 
Branching events like 10 — > 11 only become relevant for times t = 0(1/ e) since their 
rate is (3 = £ -C 1. Such processes increase the energy H = Yl,i n i an d thus require 
thermal activation. The lifetime of the resulting activated state is 0(1) because the 
excitations still coagulate with rate 7 ~ 1, either back into the original state 10 or into 
01, each with equal probability. Due to this separation of time scales the probability 
to find an activated state at any given moment t is 0(e). At low temperatures £ <C 1 
branching events therefore effectively reduce to diffusion steps 01 ^=f 10, the rate of 
which is d — [3/2 — e/2. The FA dynamics on the 0(1 /e) time scale is effectively that 
of diffusion and coagulation [32} [33l El] , 

/ \ 

—d d 7 



Wf 



d —d 7 
\ -2 7 J 



(20) 



This statement becomes exact in the limit e — > at fixed scaled time t = et. In other 
words, the evolution operators satisfy exp(W FA t/e) ~ exp(W c t/e). In units of scaled 
time t the effective master operator W c = e^W has diffusion rate d — 1/2 while the 
coagulation rate 7 ~ 1/e — > 00 becomes formally infinite. (The process nevertheless 
remains well-defined because it is diffusion limited.) 

At small but finite £ < 1 a third dynamical regime appears for times t = 0(1/ e 2 ). 
Double-branching events 100 — > 110 — > 111, which have rate 0(e 2 ), can then no longer 
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be neglected. Such triples relax with probability 1/2 into 101, so that two stable mobility 
excitations emerge. This process continuously generates new mobility excitations and 
thus counteracts the diffusion and coagulation dynamics discussed above. On a time 
scale of 0(1/ e 2 ) or larger a crossover must therefore set in that eventually drives the 
FA model into equilibrium. A simple argument shows [27] that equilibrium is attained 
only for times of 0(1 /e 3 ). 

In our analysis of aging in the FA model we focus on the nonequilibrium regime 
of times 1 <C t <C 1/e 2 or more precisely scaled times t = et = 0(1). In units of t 
the evolution is then governed by the diffusion and coagulation process Eq. (|20|) with 
rates d — 1/2 and 7 ~ 1/e 1. We will concentrate in particular on the scaling that 
emerges for large t. At finite temperature e > this scaling refers to times t ^> 1 but 
still t< 1/e since Eq. ( )20|) does not account for the onset of branching processes at 
t = 0(1/ e 2 ). In other words, whenever we talk about asymptotic behaviour for t — > 00 
it is understood that the zero-temperature limit e — > is taken first. 

We will later want to analyze response functions resulting from local perturbations 
57i = —Yli^i n i- The impact of the latter on the effective dynamics Eq. (120]) is 
easily derived. In the Glauber rates Eq. ([6]), which are defined through Wi(n) = 
[1 + exp(Aift/T)] with A { H the change in energy caused by the flip, e is then replaced 
by the local equilibrium excitation density 

£ t = 1+ Ji-H t) ,T * i < 1 + h i/n (21) 

to linear order in hi. Consequently the branching rates from sites i — 1 and i + 1 to 
site i in the FA master operator Eq. ffT9|) become (3 = ef, corresponding changes in 
the coagulation rates 7 = 1 + 0(e) are negligible at low temperature. In the effective 
dynamics Eq. ( f20|) on the 0(1/ e) time scale we then have d = f3 /2 = Ei/2 for diffusion 
from sites i — 1 and i + 1 to site % or, in scaled time t, d ~ l/2 + /tj/(2T) and 7 = 0(1/ e). 
In a linear response calculation on the t = 0(1) time scale and at low temperature e< 1, 
a field hi thus enhances the rates for diffusion towards site i. Physically, this effect arises 
because diffusion in the FA model is an activated process. 



2. A Useful Mapping 

The effective DLPC dynamics of the FA model on the 0(1/ e) time scale can be mapped 
onto a DLPA process. The latter admits an exact solution only for a special set of rates. 
We argue below that the scaling of certain observables, which includes in particular the 
domain size distribution, is robust under this modification of the rates. This will allow 
us to obtain various exact scaling results for the low temperature dynamics. We always 
focus on the 0(1/ e) time scale; to simplify the notation we will drop all tildes from now 
on, effectively using 1/e as the unit of time. 

The equivalence of entire classes of diffusion limited reaction processes has been 
established in [32] • The effective DLPC dynamics, Eq. (1201) . of the FA model in particular 
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is equivalent to the process 

/ a \ 



-d dO 
d -d 
\0 -a J 



(22) 



M+i 

which describes diffusion limited pair annihilation. Here d is the rate for diffusion 
10 01, while a is the rate for pair annihilation 11 — > 00. This result is established 
through a local stochastic similarity transformation with b = |0)(0| +2|1)(1| — |0)(1|, or 
in matrix representation [13] , 




with ft" 1 = 2 (23) 



In terms of 



one has 



3 = 6®" and £T X = (ft" 1 )®", (24) 



W c = BW^B' 1 , (25) 

if iy c and W a share the same diffusion rate d but a = 2 / y. This is readily proved: since 
the transformation is local and because W c,& = Y2i W^' 3, are pair reaction processes it is 
sufficient to consider the local master operators Wf and Wf in the jrij, rii+i) subspace. 
A simple matrix multiplication then shows that = (b ® 6)H /7 i a (6 _1 ® 6 _1 ). 

We note as an aside that there is in fact a direct mapping [34] that establishes the 
equivalence of the FA model and the DLPA process Eq. (1221) on the 0(1/ e) time scale, 
without having to proceed via the effective DLPC dynamics. The FA master operator 
Eq. ( Tl9l) exactly maps onto DLPA, Eq. ( 1221) . but with an additional pair creation process 
00 — > 11 of rate 0(e 2 ). The latter is irrelevant for the dynamics on the 0(1/ e) time 
scale so that we are once again left with the DLPA process Eq. ( f22l ). 

Let us now consider how observables transform under the mapping B. The 
dynamics of a generic observable A in DLPC is given by 

(A(t)) c = (e\Ae wct \P(0)). (26) 

By inserting identity operators BB~ l and via Eq. (1251) this becomes 

(A(t)) c = (e\BAe WH B- l \P(Q)), (27) 

where A = B' l AB. Note that (e\B = [((0| + (1|) b]® N = (e| according to Eqs. $T2]) and 
( f23l) ; likewise (e|-B _1 = (e|. In order to eliminate B~ x in Eq. (127}) an initial state |-P(0)) 
must be specified. Consider the example of an uncorrelated and homogeneous initial 
state with particle density < p < 1, 

|P(0)) = |p) = [(l-p)|0)+p|l)r JV . (28) 

Then B-^p) = [(1 - p) b^\0) + pb^\l)f N = [(1 - |p)|0) + y\l)] m = \\p) and 

(A(t)) c = (e\Ae WH \p) = (e\Ae WH \\p) = (A(t)) a . (29) 
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This equation states that the dynamics of an observable A in DLPC with initial density 
p is identical to that of A in DLPA for initial density p/2. One has to bear in mind, 
however, that A is generally a non-diagonal operator and so cannot be interpreted as a 
standard physical observable. Setting A = h, for instance, yields 

n = b- l hb= |0)(1| + |1)(1|. (30) 

The usefulness of Eq. (1291) is preserved only by the fact that A is projected onto (e|. 
For the number operator ((0| + (1|) ^ = 2(1| = 2((0| + (1|)^ as one easily verifies. More 
generally this implies [33] 

(e\n h n i2 . . . h im e WH \p) = 2 m (e\h il h i2 . . .n im e WH \\p), (31) 

for ordered (and hence distinct) %\ < %i < . . . < i m . This relates equal-time particle 
correlations of arbitrary order in DLPC and DLPA in a simple manner. A more 
convenient set of observables for our purposes will be (ii < ^) 

l 2 — 1 12 — 1 

E ilM = n ( x - "*) and = n ^ - 2hk ^ ( 32 ) 

k=i± k=ii 

which we refer to as empty and parity interval operators, respectively [5U1 EI]. This is 
for obvious reasons: E ilti2 \P(t)) = \P(t)) if the interval {h,ii + 1, . . . , ii — 1} contains 
no particles, i.e. is empty, and E il)i2 \P(t)) = otherwise. Similarly P ilti2 \P(t)) = 
(— l) m |P(t)) measures the parity of the number of particles m contained in the same 
interval. From Eq. (|29|) one has the transformation law 

( e \Eh,i2 ■ ■ ■ Ehk-iMk eW \p) = ( e l-^*i,*2 • • • Piik-ifak eW I \p)i (23) 
again assuming the %\ < 12 < ■ ■ ■ < 12k are ordered. This expression is of course 
equivalent to Eq. ( 13TT) . But its physical meaning is more appealing: empty interval 
probabilities in DLPC directly map onto parity probabilities in DLPA if the initial state 
is adjusted appropriately. 

The mapping establishes that equivalent to the effective DLPC dynamics of the FA 
model, Eq. (|2Dj) . with rates d = 1/2 and 7 ~ 1/e > 1, is the DLPA process of Eq. ([22]) 
with rates d = 1/2 and a = 27~2/e^>l. However, no exact solutions are known for 
either of these two processes. In DLPC the only exception is the process with identical 
diffusion and coagulation rates, 7 = d, for which the master operator W c admits a free 
fermion mapping [32]; obviously the same is true in the equivalent DLPA process with 
a = 2d. In the following we analyze the dynamics of these particular free fermion DLPC 
and DLPA processes in more detail. Thereafter we discuss to what extent the results 
depend on reaction rates. 

Exact results for reaction-diffusion problems are often related to integrable quantum 
spin chains [40J. Indeed, the most efficient way to generate explicit results for our free 
fermion DLPC and DLPA processes is to note that these are in turn equivalent to zero 
temperature dynamics of the Glauber Ising spin chain [52]. The latter is defined on a 
one-dimensional lattice i = 1,2,..., A of Ising spins o~i G {— 1,+1}. The probability 
p(cr, t) for being in state er = (01, 02, . . . , cr^v) at time t is governed by a master equation 
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of the form Eq. (J3J) if one replaces n by cr. Likewise the transition rates are of the form 
Eq. ([5]) but without kinetic constraint, i.e. effectively /j(er) = 1. The spin flip rates in 
the Glauber Ising model at T = are Wi(cr) = |[1 — |<7i(cXi_i + o~i + i)} and represent 
nearest neighbour ferromagnetic interactions. A connection to DLPA is established 
by considering the defect or domain wall observables rij = |(1 — <7j(x i+ i) G {0,1}. 
Writing j and j for spin configurations o = +1 and a = —1, respectively, transitions 
ttl ^ Til then correspond to diffusion 01 10 of the domain wall. The rate for this 
is d = Wi(er) = |. Similarly TIT ^ TTT translates into annihilation 11 — > 00 of domain 
walls and has rate a = Wi(cr) = 1 (which satisfies the free fermion condition a = 2d). 
The reverse process of pair creation 00 —>■ 11 is not possible under zero temperature 
dynamics. Together, these arguments show that the spin dynamics in the Glauber 
Ising chain correspond to a DLPA process of the domain walls [53] . We now exploit 
this connection to extract results for DLPA. It is convenient to introduce an operator 
formalism analogous to Eqs. ([8]) - (fill) by simply replacing the |0), |1) basis with the 
local spin basis | f), | I)- Further denote by a z the spin operator c" 2 = |T)(Tl — 1 1) (I I 
and by W a = £V (o~f — 1) Wi the master operator of the Glauber Ising chain; here 
a i = (I T)(l I + I l)(T \)i plays the role of the flip operator Eq. ( fl6l) . Using that the 
domain wall operator in the spin basis is hi = | (l — 0fer? +1 ) the parity operator for the 
number of domain walls becomes 

12—1 12— i 

*u = n - 2 "*) = n = «, (34) 

k=i\ k=i± 

since (c z ) 2 = 1. This equation expresses the obvious fact that if a ix and a i2 are aligned 
there must be an even number of sign-changes, i.e., domain walls, in the region between 
them. The fact that domain walls perform DLPA combined with Eq. (1511) supplies us 
with the result 

{e\P h , i2 ...P i2k _ uhk e WH \l P ) = (e|«...af 2fc e H/CTt |P CT (0)), (35) 

for i\ < %2 < ... < %2k as usual. The dynamics of parity probabilities in the 
DLPA process with diffusion rate d = 1/2 and annihilation rate a = 1 are identical 
to spin correlations in the zero temperature dynamics of the Glauber Ising chain. 
The initial state of the spin system IP^O)) corresponding to the uncorrelated and 
homogeneous initial state \ \p) in DLPA is determined by setting t = in Eq. fl35l) . Then 
{e\P ilM . . . P i2k _ lMk \\p) = (e\P hM \\p) . . . (e\P i2k _ lMk \\p) and (e\P itj \\p) = (1 - p)^ 
factorizes completely since \\p) is uncorrelated and thus[|| 

k 

(e|« • • • <J P.(0)> = I^ 1 - PY 2 ^- 1 - (36) 

A=l 

| Here the thermodynamic limit N — * oo was taken. In a finite system DLPA conserves the parity 
Pi, at °f the total particle number. The sectors of even and odd parity evolve independently and map 
onto Ising systems with periodic and antiperiodic boundary conditions, respectively. In the N — > oo 
limit the two sets of boundary conditions give identical results and so do not have to be distinguished. 
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These are precisely the equilibrium correlations, |P CT (0)) = |P eq (T )), for the Ising 
Hamiltonian 7i = —JY2i CT i 'i+i a ^ temperature T if 

tanh( J/Tq) = 1 - p. (37) 

Altogether, from Eqs. (133}) and (|35j) . empty interval probabilities in the DLPC process 
with identical diffusion d = 1/2 and coagulation 7 = 1/2 rates are thus given by 

{e\E iui2 . . . E i2k _ 1>i2k e WH \p) = (e|« . . . < fe e^|P eq (T )), (38) 

where i\ < i 2 < ■ ■ ■ < i 2 k- Equation (138]) is the central result of this section. It directly 
relates spin correlations in the Glauber Ising chain after a quench from temperature 
To > to T = to the dynamics of (multi) empty interval observables in DLPC for 
initial density < p < 1. The required spin correlations appearing on the right hand 
side of Eq. (138]) were derived in [54] ■ The simplest result is obtained for the initially 
filled state p = 1 which corresponds to T = oo according to Eq. fl37|) . One then has [54] 

k 

(e\E iui2 ...E i2k _ ui2k e wct \l)= £ (-1)'^^^^). (39) 

neV(k) A=l 

Here tt : {1,2,..., 2k} \— > {1,2,..., 2fc} denotes an index permutation, and (— l) 7r its 
sign; the sum runs over all permutations in the set V(k) of ordered pairings of the 
numbers {1,2, ... ,2k}. In an ordered pairing the entries are ordered within pairs, 
7r(l) < 7r(2), 7t(3) < 7r(4), . . . , 7i(2k — 1) < 7r(2fc), and the pairs are arranged in order of 
their first entry so that 7r(l) < 7r(3) < . . . < i\{2k — 1). An expression for the function 
H n (t) is given in Eq. f]A.8j) . We remark that Eq. (139]) also applies in the case p < 1, i.e. 
T < oo, if one replaces H n by the function H" given in [Si] . In subsequent sections we 
will repeatedly make use of the result Eq. (]39j) for the special cases fc = 1, 2 where it 
reduces to 

(e\E lut2 e WH \l) = H l2 . ll (2t), (40) 
^6|-E , j 1> j 2 -£/j 3j j 4 e |1) [-/T 2 — ii Hi 4 — j 3 Hi^—i^Hi^—^ -\- Hi 4 —i 1 Hi^—i 2 ^ (2t). (^1) 

Both expressions, like Eq. (139]) . only apply for strictly ordered indices. In Eq. ( ]4T]) we 
have introduced the short hand notation [ • ] (x) to indicate that all functions enclosed 
in the square brackets have the same argument x. 

2.1. Domain Size Distribution 

At this point we have the exact result Eq. ( 139]) for empty interval densities in the free 
fermion DLPC process, with rates d = 7 = 1/2. However, we are actually interested 
in the effective FA process which has the correct diffusion rate d — 1/2 but a much 
larger coagulation rate 7 ~ 1/e ^> 1. In order to understand how the latter affects 
the dynamics of DLPC it is instructive to consider the domain size distribution. Using 
spatial homogeneity we define the density of domains of size k = j — i > 1 as 

D^it) = (e\ni(l - n l+1 ) •■■(!- h^h, e WH \l). (42) 
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Equivalently Dj^t) = (e| {E i+ \j — E^j — E i+ ij + i + E it j +1 ) e^^l) as follows by re- 
writing hi = 1 — (1 — fii) and similarly for hj. In the free fermion case these empty 
interval densities are given by Eq. f HUl) so that D k {t) = [Hk-i — 2H k + H k+ f\{2t). 
Applying the recursion Eq. (1A.9j) turns this expression into 

D k (t)=e- 2t [I k ^-I k+1 ](2t) } (43) 

where the I n (t) are modified Bessel functions, see Eq. (IA.2I) . This remarkably simple 
exact result is well known [55J. Since one particle may be associated with each domain 
the exact particle concentration c(t) = (ni{t)) follows by summing this over k. The sum 
is telescopic and hence 



oo 

c(t) = J2D k (t)=e- 2t [I + h}(2t). (44) 

k=l 

Obviously the results Eqs. (1431) and ( l44l) apply exactly only for the free fermion 
(d = 7 = 1/2) DLPC process with initial density p = 1. However, it is known that 
both the concentration and the domain size distribution have remarkably robust scalings 
[56] . From Eq. ( 1441) and the asymptotic expansion Eq. ( 1A.6I) the concentration scales 
like c(t) ~ \/\fwi for t — > oo. This leading order result applies in fact to any DLPC 
process with d — 1/2 and arbitrary 7 > and regardless of the initial density p > 
|56j . To understand why, consider the scaling of the domain size distribution Eq. fj4*3|) . 
In the limit t — > 00 and for domain sizes k = 0(\fi) one finds, via Eqs. ( 1A.4I) and ( 1A.7I) . 
the scaling 

D k (t) = k -e-*I k (2t) ~ -±- ±=e-^\ (45) 
t ^/nt 2y/t 

Typical domains have a size of k = O (Vi) and density O (t^ 1 ). Integrating Eq. ()45l) over 
k = kj (2\ft) G ]0, oo[ shows that only typical domains of length k = O (Vt) contribute 
to the scaling c(t) ~ l/y/nt. Large domains k ^> Vt are exponentially unlikely and 
also small domains k \fi are of sub dominant density. In fact, the density of finite 
domains with k = Oil) is D k {t) ~ k/(2^t 3 ^ 2 ) = 0(t~ 3 / 2 ), while for typical domains 
D k (t) = (9(t _1 ). The scaling of DLPC is thus controlled by typical domains. For 
such domains to disappear, the delimiting particles first have to diffuse for a time 0(t) 
before there is a significant chance that they will occupy adjacent sites. The coagulation 
reaction itself then occurs on an 0(1) time scale for any 7 > 0. It therefore only has a 
subdominant effect on the rate at which typical domains disappear, and the process is 
said to be in the diffusion controlled regime |56j . The scaling for the density of typical 
domains Eq. (|45|) thus applies to any DLPC process with d = 1/2 and 7 > 0, and 
regardless of the initial density p > 0. The precise value of the coagulation rate 7 only 
shows up in the coefficient for the density of small domains k = 0(1) which, in any 
case, is 0(t~ 3 / 2 ). 

The above discussion illustrates that the initial density p as well as the coagulation 
rate 7 only have a subdominant effect on the evolution of DLPC [56] for t ^> 1; recall 
that in the FA model this long-time regime corresponds to unsealed times much larger 
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than 1/e but much smaller 1/e 2 . The long-time behaviour of the effective FA process 
Eq. (120]) can therefore be extracted from the analytically tractable free fermion case 
(d = 7 = 1/2). We also exploit the irrelevance of the initial density for the leading 
order long-time results to set p = 1; this simplifies the mathematics and produces 
differences to the case p < 1 only in subdominant terms. 

Note that the above simplifications will hold for all observables except for those 
which show particular sensitivity to the density of small domains; we will return to this 
point in our discussion of response functions. 

3. Two-Time Quantities 

In the following we derive exact scaling expressions for the two-time correlation and 
response functions of the FA model in the long-time regime discussed above, by 
exploiting the equivalence to the free fermion DLPC process with initial density p — 1. 

The mapping from DLPC to DLPA, which is our key for deriving results, cannot 
be applied in a practical way to multi-time quantities. Consider, for instance, 

(A 1 (t)A 2 (t w )) c = (e|Ae WCA %e^|p), (46) 

where t w is the waiting time since preparation of the system and At = t — t w is the 
time interval to the later measurement at time t. Inserting identity operators BB^ 1 
yields (e\Ai e w&M A 2 e WHv \\p) ■ But as discussed below Eq. ( 1291) the mapped observables 
A\$ — B~ l Ai 2 B are generally non-diagonal operators. While the projection (e\Ai may 
reduce A\ to a physical observable this is no longer the case for A 2 which is sandwiched 
between evolution operators. 

We now demonstrate how the mapping can nevertheless be applied to multi-time 
quantities and in a physically sensible manner if the evolution operator is known. The 
discussion will proceed in the reverse order of Sec. [2J we first consider spin dynamics 
in the Glauber Ising chain, relate these to parity intervals in DLPA and then map onto 
DLPC. We exploit the following general result for the evolution of two-spin correlations 
Fi 2 \t) = ( e \ a i 1 a i 2 eW ' 7t \' l l J } i n the Glauber Ising model at zero temperature: in the 
notation of [54j, 

if >(f) = J2 G 8(*)4 2) + H i2 - h (2t)AW, (47) 

jl<32 

where 

Gg(t) = e" 2 ' " In-nh^n] (t). (48) 

and the Aj = ( e \°~ Z j 1 °~ Z j 2 \' l l ) ) are the correlations at time while A^> = (e\ip). Eq. (T47j) 
expresses a two-spin correlation at time t in terms of those in the general state at 
t = 0. The corresponding two-spin Green's function G\j(t) is expressed in Eq. (1481) 
in terms of modified Bessel functions I n (t). The definition of the latter is recalled in 
Eq. (IA.2I) . while H n (t) is stated explicitly in Eq. flA.8j) . Equation (T47j) applies for i\ < i 2 
and it is our convention throughout that sums as in Eq. fT4T|) are taken over ordered 
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integer pairs j = (ji,j2)- Now, given that Eq. (14TI) applies for arbitrary initial states 
we may equivalently write 

<e|«e^ = G>g(t)< e K4 + H l2 ^(2t)(e\. (49) 

We refer to this object as the two-spin propagator; more general multi-spin propagators 
follow likewise from the general result given in [51] . Equation (09]) has a direct analogue 
in DLPA. Recall that according to Eq. (I341) the operator of of = measures the 

parity of the number of domain walls between sites i\ and i<i- This, combined with the 
fact that Glauber dynamics of the spin system correspond to DLPA of the domain walls, 
implies 

{e\P i e WH = G%(t)(e\P j + H l2 ^(2t)(el (50) 
h<j2 

which is the parity propagator in DLPA. A mapping to DLPC is now straightforward. 
Multiplying by B' 1 from the right and inserting identity operators B~ X B immediately 
yields the empty interval propagator in DLPC, 

(e\E ie wct =Y Gg(t)(e|^ + ^ 2 _ n (2t)(e|. (51) 

The expressions Eqs. (T49H5T]) are rather remarkable. They show that propagating a two- 
spin operator forward in time in the Glauber Ising chain is equivalent to propagating a 
parity operator in DLPA or an empty interval operator in DLPC. Nevertheless, as we 
will see, the two-time quantities that follow from Eqs. (150]) and (!5"Ii) are not related in 
a simple way. Note also that in contrast to Eq. (146j) no non-diagonal operators arise 
explicitly when the mapping is written in terms of the propagators: since Eqs. (150]) and 
( |5T]) are explicit representations of the evolution operator, a projection state (e| is always 
available to absorb contributions from off-diagonal components of mapped operators. 
We finally add that via the identity Eq. (1A.12I) each of the above propagators can be 
reexpressed in a somewhat different but equivalent form, for instance 

(e\E ie wct = (e\ + Gf](t) (e\ {E j - 1) . (52) 
3.1. Two- Time Correlations 

Using the DLPC propagator derived above we can now analyse the connected two-time 
correlation functions of mobility excitations (or defects) in the FA model, 

^•(Mw) = (el^e^^-e^ll) - (e|n i e vyc '|l)(e|%e wctw |l). (53) 

To understand nonequilibrium fluctuation-dissipation relations one also requires the 
associated response functions; these will be considered in the next subsection. It is 
useful to introduce the more general correlator 

Cij(t,t w ) = (el^e^^-e^ll) - <e|^ e^ c * 1 1> (e| e wc ^ 1 1> , (54) 
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which reduces to Cj_j(t, t w ) for i = + 1) and j = + 1) since = 1 — E ii+ i. 
The first step in evaluating Eq. (l54j) consists in decomposing the evolution operator 
e w c t _ e w c At e w c tw anc j substituting the empty interval propagator Eq. ( 1521) for 
both occurrences of (e\Eie WCAt in Eq. (|5^|) . Simplifications occur due to probability 
conservation, (e\e w ° = (e|, and the normalization (e|l) = 1 of the initial state, giving 

Cij(t,t w ) = G fl( At ) {(e\E k Eje wct -\l) - (e|£ fc e^l>(e|^-e^|l>) . (55) 

kl<k 2 

The connected two-time empty interval correlations are now reduced to connected one- 
time averages. We remark that precisely the same expression applies in DLPA if one 
replaces the empty by parity interval operators. Nevertheless, as indicated above, the 
resulting two-time correlations are different. This is for the following reason: the sum 
in Eq. ( 155]) runs over all ordered index pairs k and thus there are inevitably overlaps of 
Ek with Ej. In such cases the empty interval operators merge since (1 — fii) 2 = 1 — hi. 
Parity interval operators, on the other hand, cancel each other in regions of overlap 
because (1 — 2hi) 2 = 1. Differences in two-time correlations between DLPC and DLPA 
therefore arise due to the different reduction properties of empty and parity interval 
operators. For the case of empty intervals that we are interested in here one has 





h < k 2 < ji < 32 


Ek!,j 2 


h < Ji <k 2 < j 2 




for k x < ji < j 2 < k 2 


< E- ■ 


jx <h<k 2 < j 2 


EjiM 


ji < h < j 2 < k 2 


k Ej\,j 2 Ek li k 2 


Ji < 32 <h<k 2 



(56) 



To proceed with the derivation of Cij(t,t w ) the summation in Eq. fl55l) over the pair- 
correlation term must be broken up according to Eq. ( 1561) . with corresponding reductions 
in each sector. For the remaining averages, which are now of the form (e\Ei lt i 2 e wctw \l) 
and {e\Ei lt i 2 Ei 3 j A e w ° tv |1) with strictly ordered indices l\ < l 2 < h < U, we substitute 
the results Eqs. (TJO]) and (jHJ) from above. At this point one has a rather bulky expression 
for Cij(t,t w ). To simplify matters we focus on j = + 1) - the case of interest - 
which narrows down several summation ranges, c.f. Eq. ( 1561) . It is then a tedious but 
trivial task to show that the various sums can be regrouped into the compact form 

Ci,(j,j+i){t,t w ) = [1 - Fi(2t w )] ^2 G i,h(^ t ) H k 2 -ki(^tw) 

ki<j k 2 >j 

+ E E G%l(At)H j+1 „ kl (2t w ) [S k2tj + H^ k2 (2t w )] 

ki<j k 2 <j 

+ E E G^l(At)H k2 ^(2t w ) [6 klJ+1 + i? fcl _ i _ 1 (2t w )] . (57) 

ki>j k 2 >j 

Here antisymmetry of G^(At) under permutation of k\ ^ k 2 and H (2t w ) = was 
used. Finally we set i = (i, % + 1) and rearrange Eq. ( 1571) as follows: in each line the 
summation variables are changed from hi, k 2 to p, q > and the indices of the Green's 
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functions are simplified based on translational G\j t {At) = G\f lk+l (At) and reflection 
Gf k (At) = G_\ _ k (At) symmetry. In terms of n = % — j the last two lines of Eq. (1571) 
thereby become 

oo 

E G Un,±n-i), M ( At ) H P+1 {2t w ) [5 q , + H 9 (2Q] , (58) 

p,q=0 

where the two signs +n and — n correspond to the third and second line of Eq. (157|) . 
respectively. For later convenience we point out that hidden in this expression lies a 
leading order cancellation in the limit of large t w . Using that the Green's function in 
Eq. fl58|) changes sign when p and q are interchanged we can write 

oo 

E G Un,±n-i),( P , q) ( At ) Vpfl + H P (2Q) [S qfl + H q (2t w )] = 0. (59) 

p,q=0 

Subtracting this vanishing expression from Eq. fl58|) replaces H p+ i(2t w ) with H p+ i(2t w ) — 
5 Pt o - H p (2t w ) = -e~ 2tw [I p + I p+1 ](2t w ), see Eq. ( 1A.9I) . Similarly, in the first line of 
Eq. (jSZD we have 1 - iJ x (2t w ) = e- 2 ' w [/ + h](2t w ). Altogether our result for the 
connected two-time mobility correlation function defined in Eq. (1531) becomes 

oo 

C n (t,t w ) =e- 2t -[/ + / 1 ](2t w ) G { ^ n+lU _ M+1) (At) H p+q+1 (2t w ) 

p,q=o 

oo 

- E G (l, + n-i),( M) (^) «- 2tw [/ P + Vi](2t w ) [S q , + H q (2t w )} 

p,q=0 

oo 

- E i),(pJ At ) e_2iw [ J P + Vi](2tw) ko + H q (2t w )] . (60) 

p,q=0 

This expression is exact for free-fermion DLPC with unit initial density, but as argued 
above its scaling behaviour at large t w applies to DLPC processes with arbitrary 
coagulation rate and initial density. This includes in particular the effective FA dynamics 
we are interested in. Equation fl60|) forms the basis for our subsequent discussion; the 
leading order cancellation with respect to t w having been eliminated, it is well suited 
for long-time expansions. 



3.2. Two-Time Response Functions 

We now turn to the two-time response function conjugate to the two-time correlation 
C n (t,t w ) discussed above, which is given by 

Ri-i&U) = (eln.e^V^-ll), (61) 

where the operator Vj = T [dW c (h)/dh] \h=o accounts for the perturbation 5H = —hhj 
of the effective dynamics at t w . Note that the response contains a factor of T consistent 
with our definition Eq. ([2]). By analogy to the correlations it is convenient to consider 
the more general response 

Ri d (ft w ) = -( e \Eie WCAt V 3 e WH -\l), (62) 
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which reduces to Eq. (IBTl) for i = + 1) as is obvious from E iji+1 = 1 — hi and 
probability conservation (e\e WCAt Vj = (e\Vj = 0. As before, substitution of the empty 
interval propagator Eq. ( 1521) reduces the two-time average in Eq. ( |62i) to a one-time 
average, which is most conveniently expressed in the form 

Ri,j(t,t w ) = G^l(At)Rkj(t w , t w ), (63) 

fcl</C2 

with the instantaneous response 

Ri d (t w ,t w ) = -(el^-e^ll). (64) 

In what follows we will first concentrate on the instantaneous response; from this the 
two-time response Ri t j(t,t w ) then follows immediately via Eq. ( |63l) . The effect of the 
perturbation encoded in Vj has the explicit form 

Vj = ~(Fj-iFj - l)Vi(l - bj) + - 1) " (65) 

Here the first term is nonzero only if site j is empty and site j — 1 occupied, reflecting 
the increased rate for diffusion from j — 1 to j as discussed at the end of Sec. the 
second term similarly captures the increased diffusion from j + 1 to j. It will become 
clear below that a decomposition of Vj into an asymmetric and symmetric contribution, 
Vj = Vj + Vj , is very useful: 

V t ] = 7(^-1^ - !)(Vi " bj) + \{Fj +l Fj - l)(n i+ i - hj), (66) 



^\ j ± j 4 
yf° = ^(i^F, - l)(Vi - %) 2 + " " %) 2 - (67) 



(a) 

The operator Vy represents increased local diffusion rates j — 1 —>■ j <— j+1 for entering 
site j but decreased ones for leaving it j — 1 <— j — ► j + 1, thus ly traps diffusing 
defects at site j. Its symmetric counter-part , on the other hand, corresponds to 
increased local diffusion rates between sites j — 1 ^ j ' ^ j + 1 in both directions and 
hence speeds up the local dynamics. Accordingly, the defect response function in the FA 
model may be decomposed into the contributions 

R n (ft w ) = R^(ft w ) + R^(ft w ). (68) 

3.2.1. Asymmetric Perturbation We first discuss the effect of the asymmetric part of 
the perturbation. From Eq. fl64|) the corresponding instantaneous response function is 

R^(t w ,t w ) = -(elEiV^e^ll). (69) 

If Vj acts entirely inside or outside of Ei, then the instantaneous response vanishes: the 
empty interval density (e|.E^e * w |l) is only affected by the perturbation when applied 
at either boundary of Ei, increasing/decreasing the rate for defects to enter/exit the 
interval. Substituting the expressions ( 1321) and (1661) for Ei and V- , respectively, and 



using the commutation relation (1 — h^Fi = Fihi one shows that 



R%(U,U) = - (5 jM - Sj A ^) (e| (h^t - h h Y E il+1A2 e w 
1 

"I 



4 

-7 (kh - kh-i) (n i2 _! - h i2 ) 2 e WH -\l). (70) 
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This makes perfect sense: the Kronecker deltas enforce that j must be a relevant site and 
also carry the sign of the response, while the averages give its magnitude. In the first 
line, for instance, only the configurations O^-il^O^+i . . . Oj 2 _i and l^-iO^O^+i . . . Oj 2 _i 
contribute, since these are susceptible to perturbations of the local diffusion rates 
between sites %\ — 1 and i\ and affect the empty interval observable Ei. To evaluate 
the averages it is convenient to rewrite (n^-i — ^ij 2 = fi^-i (1 — n^) + — n^-ih^ 
whereby 

( n i 1 -i ~ n h) Ei 1+ i t i 2 = n^—iE^fa + rii x Ei^+\^ 2 — rii x -\rii x Ei^+\^. (71) 

The first two terms on the right hand side of this equation are related to the density 
of domains: the observable n^E^+i^, for instance, measures the density of domains 
li 1 0i 1 +i . . . Oj 1+n _ilj 1+n of any size n > %2 — h and hence 

CO 

(e\h h E il+1>i2 e WH -\l) = J2 D h-h+k(U)- (72) 

k=0 

Altogether we may rewrite Eq. (ITUl) in the form 

1 



(^ w ' ^ w ) 4.^*^ LJ i-2. — M V \v ^ " — 

k=l 



+ Ag(t w ,t w ), (73) 



where Aij = 5 jtil - t^-i - 5 j)i2 + t^-i and 

A^j(t w ,t w ) = ~ (S jih - ^n-i) (eln^-in^^+i^ e WHv \l) 

+\ ihh ~ Sjfr-i) (elEi^^n^h^ e WH -\l). (74) 

The function A ( ^(t w ,t w ) in Eq. ^ was isolated since it is subdominant at large t w 
and thus unimportant for the long-time scaling. The remaining terms on the r.h.s. 
of Eq. (173"]) contain a sum over the density of domains D n (t w ) which picks up its 
main contributions from domains of size n = 0{\/t^). In this regime D n (t w ) has a 
robust scaling form Eq. (j4"5j) that is independent of the coagulation rate and the same 
therefore applies to the leading long-time behaviour of the instantaneous asymmetric 
response. It is clear from Eq. (!7T|) that A^(t w ,t w ) corrects for the fact that, e.g., 
the configuration lj 1 _ilj 1 0j 1 +i . . . Oj 2 _i does not respond to perturbations of the local 
diffusion rates between sites i\ — 1 and i\. At long times such states with neighbouring 
defects have subdominant density in DLPC compared to states with isolated defects, 
and this is the intuitive reason why A^ a j(t w , t w ), whose scaling does depend on the 
precise coagulation rate, is subdominant compared to the leading term in Eq. (!73|) . (In 
the effective FA dynamics the probability of states containing adjacent defects is further 
suppressed by the fact that coagulation is essentially instantaneous.) Further discussion 



of A^(t w ,t w ) may be found in |Appendix B[ here we focus on Eq. fl73|) which controls 



the asymptotic long-time scaling. Substituting the DLPC result Eq. (I4"3l) for D n (t w ) 
produces a telescopic sum in Eq. ([751) and gives 

^ij(^w)^w) = ^Ai^e 2 * w [-^ia— *i— i + + -^ 2 -ii+i](2iw) + A^J(t w ,t w ). (75) 
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The associated two-time response function is readily obtained via Eq. (I6"5j) . The ordered 
double sum over k\ < k 2 reduces to semi-infinite ones due to Aij but the latter may be 
recombined into full summations over all integers and evaluated using the convolution 
property Eq. (1A.5j) . The resulting expression for the asymmetric part of the two-time 
response function Eq. ( 16TT) with n = i — j is most compactly expressed as 



#ij(Ww) = j (Sj,n + ^,u-i) (e\ (n n -i ~ O Eh+i,h e^ll) 



R { :\t,t w ) = d t J^~ 2t I n (t - t w )[/ n _ 1 + 2I n + I n+1 ](t + t w ) + A$(t,t w ). (76) 

3.2.2. Symmetric Perturbation The derivation of the response to a symmetric 
perturbation is largely analogous to the asymmetric case discussed above. We define 

«S(*w.*w) = -(ei^ (s) e^|l). (77) 

Again a non-zero response occurs only if the perturbation is applied at either boundary 
of the empty interval Ei. Specifically one finds 
1 
4 

+\ ihn + hn-i) (e|#M 2 -i (K-i - K) e^ Ciw |l), (78) 

which is the symmetric counterpart to Eq. (17D1 . This time we rewrite n^-i — n il = 
n^^i (1 — n^) — + n^-in^ so that 

(^u-i — ^ii) E il+ i t i 2 = n ix _\Ei x ^ — h il E il+ i^ 2 + n il _ 1 nj 1 i?j 1+li j 2 , (79) 

and once more use Eq. (1721) to express the averages (e|nj 1 £ , il+ i )i2 e^^ll) etc. in terms 
of domain densities D n (t w ). One finds 

1 

4" 

where B id = Sj^ + S jtil -i + S j>i2 + 5j,i 2 -i and 



Rfj{tw,t w ) — — --BijAj-j^tw) + A^.(i W) i w ), (80) 



A^(t w ,t w ) = +- {8 j>h + Sj^) (e|n il _in il £ il+ i, ia e WH -\l) 

+\ ( 6 j,i2 + <Wi) (e|^ ilii2 _in i2 _in i2 e WHv, \l). (81) 

Due to the difference in signs between Eqs. ( ITTj) and ( 1791) we end up with a single domain 
density D^^it^) in the symmetric response Eq. flHUj) rather than the sum we had in the 
asymmetric case Eq. ( J73|) . This has important consequences for the scaling behaviour 
of R { ^(t w ,t w ). The excess term Ag(* w ,t w ), Eq. (JSU), differs from Ag(f w ,t w ), Eq. flTJ, 
only in overall signs; we therefore discuss both quantities together in Appendix B At 
fixed i 2 — h and large t w we have, for instance, that both A^ a j(t w , t w ) and A^ s j(t w , t w ) scale 
like (9 (l/t w 3 / 2 ). This is a subdominant contribution to R^j(t w , t w ) = O (l/tw 1 / 2 ), while 
in Rfj(t w , t w ) = O (l/t w 3 / 2 ) from Eq. (IHOl it contributes to leading order. Therefore the 
scaling of i^ 8 j(t w , t w ) for small «i > 1 and large t w depends on the precise coagulation 
rate and so is different in the effective FA model and e.g. free fermion DLPC Only in 
the regime i2~i\ 3> 1 does Rf^(t W} t w ) assume a robust scaling form that is independent 
of the underlying coagulation rate. This is because the density of domains D n (t w ) has 
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a robust scaling in this regime while at the same time A^j(t w , t w ) becomes negligible in 



comparison (see Appendix B ) . Having clarified the range of validity of Eq. ( IHUl) we now 



substitute the free fermion DLPC result Eq. (T4"3"I) for the domain density, 



#g(Ww) = —B^e-^if,-^ - J ia _ il+a ](2t w ) + Ag(t w ,t w ). (82) 



Inserting this expression into Eq. ( 1631) produces the associated two-time response 
function. As for the asymmetric case the /c-summations may be carried out using the 
convolution property Eq. (1A.5|) of the modified Bessel functions. The symmetric part 
of the two-time response function Eq. (1611) with n = i — j becomes 



R^\t,t w ) = --e~ 2t {I n {t - t w )[-I n _ 2 + 2/ n - I n+2 ]{t + t w ) 

+ [/ ft _i - I n+1 ](t - t w )[/ n _i - I n+ i](t + t w )} + A^(t,t w ). (83) 

This result is exact for free fermion DLPC; its leading order scaling at large t w applies 
more generally for arbitrary coagulation rates as long as At ^> 1. The latter holds true 
since for At 3> 1 the Green's function in Eq. ( 163|) picks up its leading contributions from 
&2 — k\ = 0(\/~A~t), and in this regime Eq. ([82"]) has a robust long-time scaling form. 



4. Local Correlation and Response 

In this section we study the nonequilibrium FDT in the FA model associated with 
the local defect observable A — rij. Before presenting results for the FDR we discuss 
the dynamics of the corresponding connected two-time local defect autocorrelation and 
response functions, denoted by C&{t,t w ) and R^(t,t w ), respectively. 



4-1. Defect Autocorrelation 

The defect autocorrelation function C^(t, t w ) = (hi(t)ni(t w )) — (hi(t)) (nj(t w )) = C (t,t w ) 
is given by Eq. (1601) evaluated for n = 0, 

oo 

C d (Mw) =e-^[I + h}(2t w ) G^ lU ^ q+1) (At) H p+q+1 (2t w ) 

p,g=0 

oo 

- 2 E G (o!-i),( P ,,)(^) e " 2tw [ J P + Vi](2tw) [5 9 ,o + ^(2t w )] . (84) 

p,q=0 

Equation fl84^) is an exact expression for the connected two-time particle autocorrelation 
in the free fermion DLPC process with filled initial state. Asymptotically, for t w — > oo, 
it applies for DLPC processes with arbitrary coagulation rate and homogeneous initial 
states with finite particle density. We add that Eq. ( |84l is then exact to leading order 
in t w for arbitrary At > 0, i.e. including small At = 0(1). This is because at large 
t w , the correlation function for At of 0(1) is to leading order unaffected by coagulation 
processes, the typical distances of 0(y/t^) between particles being too large for them 
to meet. For At = 0(t w ), on the other hand, coagulation events do play a role but 
the precise extent of their 0(1) durations is irrelevant. Either way, the dynamics of the 
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DLPC process is essentially controlled by diffusion and independent of the coagulation 
rate [56J. 

Let us now discuss some limiting cases of Eq. (l8ll) . For At = 0, one has 
^i 2 j(0) = ^ii-ji^ia-ja — ^1-32^2-31 an d so om y the P = <7 = term from the first sum 
remains while the second one vanishes altogether. Using H\(2t w ) = 1— e~ 2 * w [/o+A] (2t w ), 
see Eq. (1A.9I) . we obtain 

C d (t w ,t w ) =c(t w ) -c 2 (t w ), (85) 

with c(t w ) = {rii(t w )) the defect concentration given in Eq. (jUJ). This result is 
consistent with the definition of C&(t, t w ) as it must be. In the quasi-equilibrium regime 
At > fixed and £ w — > 00, convergence of the sums in Eq. (J5"4l) is guaranteed by 
the Green's functions, with relevant contributions originating from p,q = 0(1). We 
may therefore substitute the asymptotic expansions Eqs. f)A.6j) and (I A. 101) for the in- 
dependent functions, 

-. CO 

Cd(Mw) ~^=E [ G (o!i),(- M+ i)(^) - 2Gg;_ 1)>M (At)l . (86) 
v p,q=o 

After substitution of Eq. fj48l) for the Green's functions the double sum factorizes. There 
are several cancellations and only one sum of the type Eq. (IA.1I) remains. We find 

C d (t,t w ) ~ ^=e- A */ (At) ~ c(t w ) Pr (At). (87) 

Here we have identified the probability of return to the origin for a standard random 
walker, p r (r) = e~ T Iq(t). The leading contribution to the two-time correlation arises 
from Cd(t, t w ) ~ (ni(t)ni(t w )) and its value is thus given by the probability c(t w ) of 
having a defect at site % at time t w multiplied by the probability p r (t — t w ) for this defect 
to return and occupy the same site at time t. This reasoning confirms the intuition 
that in the quasi-equilibrium regime defects are too far from each other to meet and 
react and so behave as isolated random walkers; this argument of course applies to both 
DLPC and DLPA [2Ql E7]. 

Finally we consider the aging regime At,t w ^> 1 with the time ratio At/t w fixed 
and finite. Here asymptotic expansions of the Green's functions in Eq. ( 1H11) are required, 

G (o!i),(- M+ i)( A *) = e " 2A * ihh - W*+d ( A *)> (88) 
Gg_ 1)i(M) (At) = e~ 2At [I p I q+1 - I p+1 I q ] (At). (89) 

Some care has to be taken since the leading contributions of the modified Bessel functions 
cancel. To obtain the first subdominant correction we rewrite 

I p+1 (At) = + J p+ i](At) - ^[Vi - / P+ i](At) = d At I p (At) - |^(At), (90) 

where in the second equality Eqs. ( 1A.3I) and ( 1A.4I) were used. The asymptotically 
dominant terms in the sums in Eq. fl84l) satisfy p 2 / At = 0(1). We then note that 
from the scaling Eq. (lA~7l) one has d At I P (At) = I p (At) [1 + 0(1/ At)] and, via Eq. J2QJ), 
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I p+1 (At) = Ip(At) [1 -p/At + 0(1/ At)}. From this the leading terms in the Green's 



functions are 



G 



(o!l),(-p,g+l)( At ) 
(At) 



p + q 

At 
P-Q 

At 



e-^I p {At)I q {At), 
e- 2At I p (At)I q (At). 



(91) 
(92) 



The aging expansion of Cd(t,t w ) follows by substituting this into Eq. (154"]) and using 
the scalings Eqs. (1A.6I) . (I A. 71) and flA. 11)) for the functions I n and H n . The summations 
in Eq. (1511) then take the form (l/t w ) S p 9 >o / (p/ V^w> ?/ V^w) an d become integrals in 
the aging limit, 

if^J dy(x + y)e-^ At ^ 2+ y 2 )^x + y) 



Cd(t, U 



71 



u 



dx 



dy (x - y) e -^/ At ^ x2+ y 2 )e- x %(y) 



(93) 



7T 3 /2 At 2 

with $(■) the complementary error function as defined in Eq. (lA.lip . The integrals in 
this expression have elementary solutions [57]. One shows that 

1 



Cd{t, t v 



1 



Tit 



+ 



t + t v 

t-U 



- 1 



t-u 



t-u 



! 1 1 \v 1 t ~\~ t w t t w 

arccos -i / arccos . (94) 



7T V 2t w V t + t w 7T ^tt w t + t v 

This is the exact scaling of Cd(t,t w ) in the aging limit At,t w — > 00 at fixed and finite 
At/t w . As emphasised above, it applies not only for free fermion DLPC but any diffusion 
coagulation process with finite coagulation rate, including in particular the FA model. 
The result Eq. ( 194"]) takes a simple form in the limiting cases 

At < t w : C d {t, t w ) « .— _ , (95) 



At > t w : C d (t,t w ) 



7r v / 2Att; 

37T - 8 t w 



3^ < 96) 
As expected, the At <C t w limit Eq. (1951) matches the expansion for At ^> 1 of the 
quasi-equilibrium result Eq. (157]) . Plots of the scaling function (194]) are shown in Fig. [1] 
below and illustrate the crossover from the scaling (1951) to ()96|) . 

It is interesting to compare the above scaling of Cd(t,t w ) in DLPC to the corres- 
ponding result Cd,A(t,t w ) for DLPA. The latter was derived in [20] and reads T 



Cd,A(t, t w ) 



2tt y/r+u 



1 



1 



(97) 



Vt + U 

Plots of this are also displayed in Fig. [1] below and demonstrate that the defect 
autocorrelation functions in DLPC and DLPA are qualitatively very similar. Quanti- 
tatively, the limit forms of Eq. (1971) are 

A,<<f " :C ** u ^Wflt' (98) 

§ Note that Cd(t, t w ) in [20] is defined in terms of the observable A = afaf +1 and thus differs from the 
present choice A — hi — ~(1 — of of +1 ) by a factor 4; the same applies for response functions. 
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At > t w : C d , A (t,t w ) « (99) 

From Eqs. (195)980 we observe that Cd,A(t,t w ) « |Cd(t,t w ) in the regime At <C t w . 
This is a direct consequence of the quasi-equilibrium behaviour Eq. ( l87l) . which applies 
for DLPC as well as DLPA, and the fact that defect concentrations are related by 
(^-i(^w))a = |(^j(^w))c, c.f. Eq. f l3T|) . But in the opposite regime At 3> t w one finds 
C d , A (t,t w ) « rC d (t,t w ) with r = 37r/(67r - 16) » 3.31 from Eqs. flSSES}- Here defects 
interact and we thus expect that the nontrivial ratio r is determined by genuine many 
body effects. These must clearly differ between DLPC and DLPA, given that - amongst 
other things - the two processes have different scaling forms for their domain size 
distributions 1551. 



4-2. Local Response Function 

According to Eq. (1681) the two-time defect autoresponse function in the FA model 
may be decomposed into two contributions R d (t,t w ) = R d (t, t w ) + R d s \t,t w ). As the 
asymmetric and symmetric perturbations have very different effects on the local defect 
dynamics we discuss the corresponding response functions separately. 



4-2.1. Asymmetric Perturbation The two-time defect autoresponse to the asymmetric 
part of the perturbation is given by Eq. ( 1761) with n = 0, 

#i a) (Mw) = d tv e- 2t I (t - t w )[I + h](t + t w ) + A^ a) (t,t w ). (100) 

The correction Aq^ (t,t w ) is subdominant at large t w and for any At > and will 
therefore be neglected in the following. Due to the explicit t w -derivative in Eq. (11001) 
the corresponding step response function can be obtained directly as 

xfMw) = f dri?i a) (t, r) ~ e~ 2t {[I + h](2t) - I (t - t w )[I + I x ]{t + t w )} . (101) 

Plots of Xd\t,t w ) are shown in Fig. [1] below. This result is interesting in various ways. 
Consider first the quasi-equilibrium behaviour at fixed At > and large t w , 

xi a) (t,t w ) ~ [1 - e- At / (At)] ~ c(t w ) [1 -p r (At)} , (102) 

y7rr w 

to leading order in t w . The form of Eq. (11021) reflects a simple mechanism: at time t w 
there is a homogeneous concentration c(t w ) of defects in the system. But application 
of the asymmetric perturbation V i , which traps diffusing defects on site i, yields an 
increase of (fii(t)} as given by Eq. ( I102j) . The step response Xd^Mw) approaches c(t w ) 
within a time interval At = 0(1), see Fig. [TJ In the aging regime, where both At and 
t w are large and comparable, one has from Eq. (11011) that x d (t,t w ) ~ c(t). The step 
response then decreases along with the density of defects in the system. This explains 
the non-monotonic shape of X d (t,t w ) shown in Fig. [TJ 

It is important to note that Xd^j^w) is dominated by the small-At behaviour of 
Rd^itjtw)- Indeed, in the quasi-equilibrium regime Eq. (11021) reproduces 

4 a) (Mw) = -d tw x d a) (Mw) ~ -c(tw) d AtPr (At), (103) 
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Figure 1. Left panel: scaling form of the two-time defect autocorrelation Cd(t, i w ), 
Eq. (f94|) . in the FA model for t w = 10 2 , 10 4 and 10 6 (full lines) and, for comparison, 
the corresponding DLPA result Eq. (|97|) at the same values of t w (dashed lines). Right 
panel: plots of the step response Xd (*> ^w) to ^ ne asymmetric part of the perturbation, 
Eq. (|101[) . for t w = 10 2 , 10 4 and 10 6 from top to bottom, respectively. 



consistent with Eq. (llOOj) ; this result is positive since p r (At) is a decreasing function 
of At and hence —dAtPr(At) > 0. But in the aging regime Xd (Mw) ~ c (^) to leading 
order as discussed above, which would naively suggest R^\t,t w ) = — <9 tw c(t) = 0. More 
precisely it means that the aging part of R^\t,t w ) only contributes subdominantly to 
(t,t w ). From the full expression Eq. ( 11011) . or equivalently Eq. ( 11001) . one deduces 
that in fact 



1 

2^ 



(104) 



.Vt + u(t-t w )^ ^/t—tZ{t + t w fi\ 

For At -C t w the first term in Eq. fjl04j) dominates, matching the large At behaviour of 
the quasi-equilibrium result Eq. (11031) . For At ^> t w , on the other hand, the terms in 
Eq. (11041) cancel to leading order in t w /At; then R^\t,t w ) m t w /(n At 3 ). 

We add that a response function analogous toi?i a) (t,t w ) was calculated for DLPA in 
[20] . Denoting the latter by Rd,A(t,t w ) one has from [20] that Rd,A(t,t w ) = ^R^ (t,t w ). 
In terms of the two-time defect autoresponse for an asymmetric perturbation there is 
therefore no distinction between DLPC and DLPA, except for the factor of | which has 
its origin in (hi(t w )) a = ±(hi(t w )) c , see Eq. (J3TJ) . 



4-2.2. Symmetric Perturbation Setting n = in Eq. ([531) gives for the autoresponse to 
the symmetric part of the perturbation, 

R d s \t,t w ) = -h~ 2t Io(t-t w )[Io-I 2 ](t + t w ) + A^(t,t w ). (105) 

As discussed in Sec. 13.2.21 the behaviour for short time differences At = 0(1) of this 
result does depend on the precise coagulation rate. So in this regime Eq. (11051) applies 
to free fermion DLPC but not the FA model. This does not constitute a problem, 
however, since Eq. f ll05[) scales as O (t w ~ 3//2 ) at large t w and fixed At. It is therefore 
subdominant against the short-time scaling of the asymmetric response Eq. f l 1 3 f) . Only 
in the aging regime does the symmetric response contribute significantly to the overall 
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response function; there it is independent of the coagulation rate as explained above 
and so does apply to the FA model. Taking the appropriate limits in Eq. (I105j) gives 

R *^~-l vr^V^ - (106) 

For At < t w this simplifies to R { ^\t,t w ) « -l/[27r V / At(2t w ) 3 / 2 ] while for At > t w 
one has R^\t,t w ) « — l/(27rAt 2 ). Note the overall negative sign of the response: to 
understand this result qualitatively recall that the symmetric perturbation increases 
the local diffusion rates between sites i — 1 ±=t i ±=p i + 1. A response can only occur 
if a defect is present on one of these sites at time t w . During the time interval where 
the perturbation is applied the dynamics of this defect are accelerated. In and by itself 
this does not have any effect on the defect concentration (rk(t)) at site i. However, the 
"extra time" for the defect manifests itself through the fact that it will now coagulate 
earlier with a neighbouring defect. The local concentration is therefore reduced, and 
this explains the negative sign of the response. Physically this behaviour derives from 
the activated character of the FA dynamics [38J. The increased local diffusion rates 
in the effective FA model represent locally reduced energy barriers for the creation of 
defects. In equilibrium this would increase the defect density; out of equilibrium the 
dominant effect is a speeding up of the relaxation of the defect density towards zero, 
leading to a reduction in the number of defects. 

The scaling of the step response function X d (t, t w ) in the aging regime is obtained 
correctly by integration of Eq. ( 11061) . the reason being that small- At contributions from 
R^\t,t w ) only amount to subdominant corrections in the integral. Thus, 

xi S) (Mw) = f dr<(t,r) ~ -±^jL, (107) 

in the aging regime At, t w 3> 1 with fixed and finite At/t w . For At -C t w this reduces to 
Xd°(*,*w) « -v / At/(2 v / 27Tt w 3 / 2 ) while for At > t w one finds xJ } (Mw) « -l/(27rAt). 
Comparison of Eq. (11071) with Eq. (110 II) shows that the overall step response in the 
aging regime Xd{t,t w ) = Xd (t,t w ) + Xd (Mw) ~ c (t) is dominated by the short- At 
behaviour of R^\t,t w ) as discussed above. The aging parts of R^\t,t w ), Eq. pi]), 
and ^ s) (t,t w ), Eq. (EI]), are both of the same order but only contribute subdominantly 
in the overall step response Xd(t, t w ). The consequences of this will become clear below. 

4.3. FDR and FD plot 

Based on the above results for connected two-time defect autocorrelation and response 
functions one immediately obtains the associated fluctuation-dissipation (FD) plots. 
We use normalized FD plots (20J, [59] [60] of x(t, t w ) = x(f , *w) /C(t, t) versus 1 - C(t, t w ) 
where C(t,t w ) = C(t,t w )/C(t,t). Normalization scales the plot so that the abscissa 
always lies in the range from zero to one. It is crucial that the plot is parameterised by 
t w with t held fixed [T9J, [59]. Only then is the slope of the plot guaranteed to be X(t, t w ) 
as one readily verifies from Eq. ([3]) and R(t,t w ) = — <9( w x(t, t w ). 
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Figure 2. Left: Normalised FD plot for the local defect observable A = hi- Data 
were obtained by numerical evaluation [58] of the full expressions {53} , (|100|) , (|105[) at 
t = lO 1 / 2 , 10 1 , 10 3 / 2 , 10 2 with t w e [0,t]. For increasing t the curves move upwards, 
converging to the equilibrium line of unit slope. Inset: The data in the top-right corner 
on a magnified scale, not containing the t = 10 1 / 2 curve. Right: Plot of the FDR 
X d (t,t w ) versus t w /t for t w = 10 1 / 2 (dotted), 10\ 10 3 / 2 , 10 2 and oo (full line). The 
curves for finite t w were again obtained by numerical evaluation of the full expressions 
while the limit curve follows from the scaling forms (|94|) . (|104j) and (I10G|) . 



The FD plot for the local defect observable in the FA model is shown in Fig. [2j It 
quickly converges to the equilibrium line Xd = 1 — Cd, consistently with observations 
made in computer simulations [2H EE] . However, a plot of the FDR X d (t, t w ) , also shown 
in Fig. El clearly demonstrates that equilibrium FDT is not satisfied. In fact Xd(t,t w ) 
steadily decreases with increasing t and actually turns negative around t w /t ^1/3. This 
is reflected in the FD plot, which bends down as Cd — > 0. The plot of X d (t,t w ) also 
indicates that a nontrivial scaling form X d (t,t w ) ~ X d (t w /t) exists in the aging limit; 
we return to this below. 

In order to understand the shape of the FD plot let us first consider the dynamics in 
the quasi-equilibrium regime of large t w but finite At > 0. Here Cd(t, t w ) ~ c(t w )p r (At) 



as given by Eq. flST]) and Xd(Mw) ~ Xd (Mw) ~ c(t w )[l - p r (At)], Eq. ([102]), is 



dominated by the asymmetric perturbation. To leading order in t w we then have for 
the normalized quantities \& = 1 _ Cd = 1 — Pr(At) and equilibrium FDT is satisfied. 
Furthermore, since Cd ~ p r {At) becomes negligibly small already in the At = 0(1) 
time sector, aging effects are compressed into the top right corner of the FD plot. 
Asymptotically only the quasi-equilibrium regime is visible in the FD plot [20l [371 EH] ■ 
Consider next the aging regime where both At and t w are large, but with a finite 
ratio At/t w . The normalized correlation then is Cd = 0(1/ \fi). The dominant scaling of 
the step response Xd ~ c(i) would give Xd — 1- This is modified by aging contributions 
from i? d a) (t,t w ) and R d s \t,t w ) so that 1 - Xd = 0(l/y/i). In the FD plot, aging effects 
are only visible in a region of size 0(l/\/t) in the top right corner of the plot; this 
regime is highlighted in the inset of Fig. [21 Effects of the asymmetric perturbation 
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alone, which are accounted for in Xd (Mw)> Eq. ( llOip . would cause the plot to turn 
horizontal with X°° = 0; this would be the analog of the DLPA result in [20]. However, 
the additional negative contribution of Xd\t,t w ), Eq. (I107p . in the overall step response 
Xd(t,t w ) = Xd (^)^w) + Xd\t'tw) m akes the plot non-monotonic with X°° < 0, which 
is consistent with simulations [38J. Note that this non-monotonicity is genuine and not 
due to incorrect parameterisation of the plot with t instead of t w [T3] . 

The above discussion shows that for defect observables an FD plot gives a rather 
unsatisfactory graphical representation of the FDT violation in the FA model. It is 
more appropriate to look instead directly at the time dependence of the FDR X d (t, t w ), 
as shown in Fig. [2j The FDR involves the response functions R^\t,t w ) and R d \t,t w ). 
In contrast to the associated step response functions, these are not dominated only 
by the quasi-equilbrium dynamics at small At and so can reveal genuine aging effects. 
Looking at the relative importance of the asymmetric and symmetric contributions, 
one sees that in the aging regime At, t w ^> 1 with At/t w = 0(1) both R d a \t, t w ), 
Eq. (PUD , and i? d s) (t,t w ), Eq. ffT06]) . are of the same order £>(l/t w 2 ). Keeping in 
mind that R d (t,t w ) > while R d \t,t w ) < and looking at the plot of X d (t,t w ) 
in Figure El we notice that R d \t,t w ) in fact becomes dominant for At ^> t w ; explicitly, 
# d s) (Mw) « -l/(27rAt 2 ) while i? d a) (t,t w ) « t w /(vrAt 3 ) = 0((t w /At)\R { d \t,t w )\). This 
causes the change of sign in Rd(t, t w ) and thus Xd(t, t w ). So, while in the At <C t w 
limit the response is dominated by the trapping effect of the asymmetric part of the 
perturbation, it is the speedup effect of the local dynamics through the symmetric 
component of the perturbation that dominates for At ^> t w . This speedup is the 
consequence of coupling of the observable to energy barriers in the activated dynamics 
of the FA model. 

By working out the scaling of d tvr Cd(t,t w ) from Eq. (1911) and using the results 
Eqs. (11041 11061) one obtains the scaling form of Xd(t, t w ). We do not state the explicit 
result since it is somewhat bulky; importantly, however, its time dependence is only via 
the ratio t w /t. A plot of Xd(t w /t) is contained in Figure EJ It crosses over from Xd = 1 
in the quasi-equilibrium limit to Ad = for t 3> t w ; the asymptotic FDR is 

A d °° = - 3?r ifi « -3.307. (108) 

O 7T — lb 

We recall that the results shown in Fig. [2] are fully exact for the free fermion DLPC 
process (with filled initial state). In the limit of large t w , where the precise coagulation 
rate and initial state density become irrelevant, the results - including the scaling form 
of A d (t w /t) and the value X™ - apply to arbitrary DLPC processes and, in particular, 
to the FA model. Finite time data were included in Figure [2] only in order to given 
an indication of the speed of convergence (for free fermion DLPC) to the long-time 
asymptotic results. 
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5. Global Correlation and Response 

Given that the FA model can be viewed as a coarsening system, and following the 
philosophy of [20], we now analyse the nonequilibrium FDT for the global observable 
A = which is nothing but the energy. Our expectations are, firstly, that we 

should find the same X°° as for the local defect observable [20, [35] and, secondly, that 
a nontrivial limit FD plot should exist in this case [20l [60l [61] . As above we first study 
separately the connected two-time energy correlation and response functions C e (t,t w ) 
and R e (t,t w ), respectively; both will be normalised by 1/N to get quantities of order 
unity. Then the resulting FDR and FD plot are discussed. 



5.1. Energy Correlation 

We obtain an expression for C c (t,t w ) by summing Eq. (1601) over n since, using 
translational invariance, C e (t,t w ) = (1/N)^ i ^[(hi(t)fij(t w )) — (ni(t))(fij(t w ))] = 
^2 n C n (t,t w ). Via the convolution property Eq. ( 1A.5I) one readily shows that 

E Gr Sn + i),(-P 1 « +? )( A *) = e ~ 2A U+« - Wd(2At), (109) 

n 

J2 G Un,±n-m P J At ) = e~ 2A V P - q -i ~ WiPAt). (110) 

n 

Thus we have from Eq. f[60|) 

oo 

C e (t,t w ) = e- 2t [I + h](2t w ) [h+ q ]{2At) H p+q+1 {2t w ) 

p,q=0 

oo 

-2 e~ 2t \h- q -i - Wi](2Ai) [I p + / P+ i](2t w ) [6g, + H q {2t w )} . (Ill) 

p,q=0 

The energy correlation function C e (t,t w ) is constant - with respect to At - in the 
quasi-equilibrium regime. This is obvious when considering that defects are typically a 
distance O (y/tZ) apart and can therefore not meet within a time interval At = 0(1). 
Thus no coagulation processes are possible and the energy remains constant. It is 
somewhat awkward to derive this directly from Eq. fillip because the double sum in the 
second line only reduces to a one-dimensional sum; the assumption At = 0(1) imposes 
merely that p — q must be of 0(1) to get a significant contribution. One can get the 
result indirectly, however, by first going to the aging regime and then considering the 
limit At t w (see below). The correlation function at At = itself, which gives the 
energy fluctuations, is much simpler to work out using J n (0) = S n> o. In this case it can 
be shown that 

C e (t w , t w ) = 3e- 2 ^[/ + /i](2t w ) - e- 4 *»[3/ + ±h + / 2 ](4t w ) (112) 

An aging expansion of Eq. fillip is obtained rather easily using Eqs. (1A.4|) . (1A.6I) . (1A.7I) 
and flA.lip . In the aging limit the double sums turn into integrals and one finds 

C e(Mw)~Z^/ dx / dy(x + y) e -^ At ^ 2 <S>(x + y) 



71 



^0 
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dx / dj/fi-nje-^^M e-^d/). (113) 



Again there are elementary solutions to these integrals [57], giving 
2 



C e (t, t w ) 

in 



1 1 Vt=T v 



\ft \Jt + t w nt 



2 

4 



7T 



3/2 



1 . / tw 2 . t w 

arcsm \ , arcsm — 



, , . (114) 

t w V t \/t 4~ t w t 

This leading order long-time scaling applies to all DLPC processes with finite coagulation 
rate in the limit t w ^> 1 and for any At > 0. Coagulation rate effects are in principle 
present for At = 0(1) but are subleading for large t w because they are driven by 
the density of small domains. For At -C t w the expression ( II 14ft approaches a constant 
value as anticipated and matches up with the large-t w expansion of the equal-time value, 
Eq. (11121) . Specifically, the limit forms of Eq. (11141) are, 

^ 2\/2 

At < t w : C c (t, t w ) « _1 , (115) 

V7T t w 

At»t w :C c (t,t w )^|^^. (116) 

Figure [3] below shows plots of the scaling function Eq. fll 14j) . illustrating the crossover 
from the plateau Eq. fll 15)) to the long time scaling Eq. (11161) . 

In analogy with the discussion in Sec. 14.11 we again compare the scaling Eq. (11141) 
of C e (t, t w ) in DLPC with the corresponding result C et A(t,t w ) in DLPA. From Ref. [20J, 



CM „ _L (_i _ _y , (117) 

bearing in mind that the result in [20J contains an extra factor 4. Plots of Eq. fl 1 1 T[) are 
also included in Fig. [3] below, and cross over between the limits 

2 /o 

At < t w : C e>A (t, t w ) w 7-^=, (118) 

2v 7r t w 

At»t w : MMw)«^^. (119) 

In the regime At -C t w we can again establish a simple relationship between the 
functions C et A(t,t w ) and C c (t,t w ). It is sufficient to consider the instantaneous energy 
fluctuations for t = t w since both functions are, to leading order, independent of 
At when At < t w . Now, by definition C^ A (t w ,t w ) = (1/N) Ylijii^ii^fijit^a - 
(ni(t w ))a(%(*w))o]- Using the mapping Eq. (EI|) we have, for i ^ j, (hi(t w )nj(t w )) a = 
\(fii(t w )hj(t w )) c . However, the diagonal contributions i — j in this sum yield 
(ni(t w )hi(t w )) a = (hi(t w )) a = ^(hi{t w )) c . Keeping this in mind one readily shows that 
C c ,A(t w ,t w ) = \ [Ce(t w ,t w ) 4- c(t w )] with c(t w ) = (nj(t w )) c ~ l/v 7 ^ as usual. This 
identity is indeed satisfied by Eqs. (II 15111 18)) . While we found C^it, t w ) ~ |Cd(t, t w ) 
for defect correlations in the regime At -C t w , the energy fluctuations are related by 
C e) A{t,tw) ~ r C e (t,t w ) with r = 1 4- l/\/2 ~ 1.71. Interestingly, in the opposite 
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regime At t w we have from Eqs. (I116II119P that C 6) a(Mw) « r C e (t, t w ) with the 
same nontrivial ratio r = 3n/(6ii — 16) ~ 3.31 as for the defect correlations. IjjJ 

5.2. Energy- Temperature Response 

By similar arguments as for the correlation one sees that the normalised energy 
response to a field acting uniformly on all sites i of the system is given by R e (t,t w ) = 
En^n(Mw)- Using the decomposition of R n (t,t w ), Eq. (1681) . we again discuss the 
asymmetric, i?i a) (t,t w ) = J2 n Rn\t,t w ), and symmetric, R { e\t,t w ) = J2 n Rn\t,t w ), 
parts separately. 

5.2.1. Asymmetric Perturbation Summing R^\t,t w ), Eq. ( 1761) . over n and using the 
convolution property Eq. (1A.5I) gives 

i^ a) (Mw) = <9t w e- 2t [/ + ii](2*) + A^(t,t w ) = 0. (120) 

The first term in Eq. (11201) is independent of t w and hence the derivative vanishes 
trivially. It can be shown that the subdominant correction Ae (t, t w ) also vanishes 
exactly, so that altogether R^\t,t w ) = 0. With hindsight this result is obvious: when 
applied to all sites, the asymmetric perturbations of the local diffusion rates cancel. We 
remark that the same is true for the analogous response function R e ,A{t, tw) in DLPA 
studied in Ref. [20]. 

5.2.2. Symmetric Perturbation The sum over R ( n\t,t w ), Eq. ([MD, follows again 
straightforwardly using the convolution property Eq. ( 1A.5I) . Here we find 

R®(t,t w ) = -e~ 2t [/ - h](2t) + Af(t,t w ). (121) 

The first term is exactly d t c(t) = — e~ 2t [Io — I<2\(2t) with c(t) as given in Eq. (JHJ). This 
result is easily explained: first, since the response is normalised by 1/N it measures 
changes in the defect concentration c{t) = (1/N) ^("■itX)); the latter being just 
the normalised version of the energy. Now consider the effect of the corresponding 
symmetric perturbation. When applied locally it increases the diffusion rates between 
sites i — 1 <=> i <=> i + 1. Thus, in the global case, all diffusion rates are increased 
and the perturbation essentially speeds up the dynamics of the system. This is 
equivalent to saying that the system gains some extra time St to evolve and therefore 
Ri s \t,t w ) = [c(t + St) — c(t)}/St = d t c(t). In the aging regime this gives 

R^(t,t w ) ~ d t c(t) ~ - 2 ^ 3/2 - (122) 

| We note that the ratio r coincides with the absolute value of the asymptotic FDR X^ 3 , Eq. (|108ll . 
For defect observables this originates from the fact that c^C^a^j tw) ~ l/(27rf 2 ) in DLPA for t ^> t w , 

see Eq. (|99l) . scales in the same way as R^(t,t w ) l/(27rf 2 ) in DLPC from Eq. ([1®- Likewise 

for the energy 9t w Ce,^(t, i w ) ~ l/(2v / 7rt 3 / 2 ) in DLPA for t > t w according to Eq. (|119p matches the 
scaling of R.i s \t,t w ) ~ — l/(2y / 7rt 3 / 2 ) in DLPC from Eq. (I122j) below. However, this seems to be a 
mere mathematical coincidence with no obvious physical implications. 
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Figure 3. Left panel: scaling form of the normalized two-time energy correlation 
C* c (t,i w ), Eq. pT4]) . in the FA model for t w = 10 2 , 10 4 and 10 6 (full lines) and, 
for comparison, the corresponding DLPA result Eq. (|1 17[) at the same values of t w 
(dashed lines). Right panel: scaling plots of the energy-temperature step response 
- Xe (t,t w ) = -Xc S) (i,i w ) obtained from Eq. (fl22|) for t w = 10 2 ,10 4 and 10 6 from top 
to bottom, respectively. 



The global symmetric perturbation represents globally reduced energy barriers in the 
activated dynamics of the FA model, thus there is accelerated relaxation and a negative 
response. The scaling xi (t, tw) ~ ~~ (t — t w )/(2-y/7rt 3 / 2 ) follows immediately from 
Eq. (11221) and is shown in Fig. [3J It crosses over from X e\t,t w ) « -At/(2^t 3 J 2 ) 
for At < t w to Xc\t,t w ) w -1/(2 v / ^A7) when At > t w . 

The above reasoning applies for At 1 where the second term in Eq. (11211) is 
subdominant. For small At, it cannot be neglected, and accounts for deviations arising 
from the fact that the perturbation does not increase all rates: that is, coagulation 
rates are left unchanged. Since coagulation processes are the only ones that can give 
an instantaneous decrease of the energy, the instantaneous response must then in fact 
vanish, R^\t w ,t w ) = 0, and this can be verified from Eq. (112ip by explicitly evaluating 
A (s) (t t 1 



5.3. FDR and FD plot 

By combining the results derived in this section one obtains the FD plot and FDR 
presented in Fig. HI Consistently with our expectations [201 EH EI] a nontrivial limit 
FD plot exists for the global energy observable in the FA model. It defines the scaling 
form X e (t,t w ) ~ X e (t w /t) of the associated FDR. Since the energy-temperature step 
response Xe is negative, the FD plot lies below the x — axis. Because also the impulse 
response R e is negative the slope of the plot, and thus the energy FDR X e (t,t w ), are 
likewise negative throughout. 

As for the local defect observable it is crucial to use the correct parameterization 
of of the FD plot, with the observation time t fixed and t w G [0,t] the running 
parameter [59], ED] . Had we kept t w fixed and parameterized the plot with t, the non- 
monotonic shape of Xe(t,t w ) = Xe (t,t w ) shown in Fig. [3] would have led to a non- 
monotonic FD plot, incorrectly suggesting a sign change in the FDR A e (t,t w ). Only 
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when parameterized with t w does the slope of the FD plot correspond to the FDR X. 
The introduction of no-field methods for the measurement of step response functions 
[62] [63] has made it feasible to generate such FD plots with relative ease in numerical 
simulations. 

Visually the limit FD plot in Fig. H] appears to be a straight line, but this in fact 
not the case. The slope becomes steeper from left to right and correspondingly the FDR 
X e retains a dependence on the ratio t w /t, c.f. Fig. HI The latter is easily derived from 
the scaling of Re\t, t w ), Eq. (I122p . and that of d tvj C c (t, t w ) from Eq. (11141) . and has the 
explicit form 



Xju) 



;i + m) 3 / 2 



2 

1 H — arcsin u 

7T 



2 

+ — 

TIU 



[\ — w) 3 / 2 arcsin ^fu 



1 + u \fu 



(123) 



where u = t w /t. The slope of the limit FD plot at the origin is given by X e (t w /t — * 1) = 
-(1 + V2) « -2.41. This means that there is no quasi-equilibrium regime, a feature 
that seems to be generic for global observables, at least in coarsening systems [20]. More 
importantly, the asymptotic energy FDR X^° is obtained by taking the opposite limit 
tw/t — > 0, and this gives 

X- = ]imX e (u) = ~ 7 r^ L - « -3-307. (124) 

u^O 6 7T — 16 

So indeed X^ = X£°: the asymptotic FDRs of local and global observables are identical 
[20] [35] [38] [60] IBT] . As for the case of local defect observables our results will be exact to 
leading order for the FA model in the long-time limit; this includes all aging expansions, 
the limit FD plots and the expressions for X e (u) and X%°. Again, finite time data are 
included in Fig. H] only to give an indication of the speed of convergence (for free fermion 
DLPC) to the asymptotic scalings. 

While interesting in itself, the robustness of X°° = X^° = X£° is also rather useful 
from a practical perspective: accurate measurement of X°° in simulations based on 
the local defect observable is quite difficult as is clear from the FD plot in Fig. [2j the 
aging properties of and Xd are disguised by quasi-equilibrium contributions. For 
the global energy observable, on the other hand, a nontrivial limit FD plot exists, see 
Fig. HI and can be probed directly in simulations. Indeed, we have recently carried out 
such simulations for the FA model and the results are in very good agreement with the 
theoretical predictions set out above |38j. 



6. Summary and Discussion 

In this paper we have analysed the nonequilibrium fluctuation-dissipation (FD) 
behaviour of the one- dimensional FA model after a quench to low temperature. The 
dynamics in this regime is essentially a diffusion-coagulation process, and the link 
from this to diffusion-annihilation and thence to the Glauber-Ising chain allowed us 

% We remark that the nonuniform convergence of X e (t, t w ) in the right panel of Fig. [4]is due to short- 
time contributions from A[ s \t, t w ) in Eq. (jl21[) which yield Ri s \t w ,t w ) — and hence X e (t w ,t w ) = 0. 
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Figure 4. Left: Normalised FD plot for the energy. The curves correspond to t = 10 1 
(dotted), 10 2 , 10 3 , 10 4 and oo (full curve) with t w G [0,i\. Finite time data is obtained 
by numerical evaluation [58] of the full expressions and (]121|1 while the limit 

curve follows from Eqs. (|114|122jl . Right: Plots of X c (t,t w ) versus t w /t for £ w = 10 1 
(dotted), 10 2 , 10 3 , 10 4 and oo (full curve) obtained from the same equations. 



to obtain exact scaling results for the correlation and response functions of defects 
(mobile regions). We were thus able to resolve a puzzle posed by simulation results [21J: 
the equilibrium behaviour suggested by numerical FD plots of local quantities is only 
apparent and hides the underlying non-equilibrium effects. The root cause of this is 
that the regime where the FDR X differs from unity becomes progressively compressed, 
as times grow, into a corner of the plot. The FDR X itself is nevertheless a nontrivial 
scaling function of t w /t. A second key observation was that local and global defect 
observables give consistent values of the asymptotic FDR X°°, reinforcing the generality 
of our findings for Ising models [20l [60], [61] . 

An unexpected result of our analysis is that X°° is, in fact, negative. This 
seems counter-intuitive, certainly if one were to try to interpret T/X°° as an effective 
temperature. However, our analysis showed that the negative value of X°° can be 
traced to the behaviour of the response functions, which are sensitive to the activated 
nature of the dynamics [TF] [22] [38]. In the global response function, for example, 
the applied perturbation corresponds effectively to an increase in temperature. In 
equilibrium this should increase the energy. However, in the out-of-equilibrium dynamics 
the increased temperature speeds up the relaxation of the energy and thus actually 
decreases its value. In agreement with this intuition, we found a negative (impulse) 
response R e (t,t w ). Except at short time differences, it is t w -independent and just given 
by the time derivative (at time t) of the unperturbed energy relaxation: the small "extra 
time" provided by the perturbation has the same effect whenever it is applied. 

As regards the wider applicability of our results, we have repeatedly emphasised 
that the two-time correlation functions we have calculated are exact for generic diffusion- 
limited pair coagulation (DLPC) processes in one dimension, in the long-time scaling 
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limit. In the regime of small time differences, At = t — t w <C t w , the behaviour is 
easily rationalized in terms of non-interacting random walkers. Rather less trivial is the 
behaviour in the opposite regime t 3> t w . Here we find scaling exponents identical to 
those for diffusion-annihilation (DLPA) processes, but quantitatively different scaling 
functions. The propagators used to derive these results are exact for the free fermion 
DLPC and DLPA processes. Higher order propagators could be obtained immediately 
from the results in [51] and would form the key to analysing multi-particle and/or 
multi-time correlations in these processes for general initial conditions. 

Our response function results are similarly exact for generic DLPC processes in the 
long-time limit. Here one has to bear in mind, however, that the way the perturbation 
by the field enters can vary between models. For example, a DLPC process can be 
realized as the q — > oo limit of the one- dimensional Potts model [55] with Hamiltonian 
H = — J2i 3si,s i+1 , Si G {1, ... , q}, and Glauber rates. At low temperatures the domain 
walls rii = 1 — S SijSi+1 have DLPC dynamics with coagulation rate 1 and diffusion rate 
from bond to + i + 2) of 1/2 + (hi + h i+ i) / (4T) , with the fields hi conjugate 

to 7ij. This is precisely the asymmetric part of the perturbation that we found in the 
FA model; the symmetric part, which in the FA case contains all the activation effects, 
is absent. Because the asymmetric part of the response behaves similarly in DLPC and 
DLPA, see comments at the end of Sec. 14.2.11 and Sec. 15.2.11 the FD behaviour is then 
also like in the annihilation case [20], with an asymptotic FDR X°° = 0. 

Our results also shed light on a recently proposed method for measuring response 
functions in spin systems [M] and deriving associated nonequilibrium FD-relations. This 
approach assumes a specific asymmetric form of the perturbation, which in the FA case 
effectively amounts to discarding the symmetric part of the response. As the discussion 
of the Potts model above shows, this has drastic effects; in particular, it would lead one 
to conclude incorrectly that the FA model has a vanishing asymptotic FDR instead of the 
nontrivial value X°° = —3tt/(6tt — 16). The method of [M] and the out-of-equilibrium 
FD-relation derived from it can therefore not be accepted as generally valid. 

It is worth noting that our result for the response function of the FA model has 
recently been quoted as being consistent with the predictions of a generalized form of 
local scale invariance [65]. This is indeed so for the global response Eq. ( 11221) but not for 
the local response function Eqs. (11041)1061) . because of the factors of t + t w that appear 
in both the symmetric and asymmetric parts. This suggests to us that the agreement in 
the global case is coincidental and cannot be taken to imply that the FA model genuinely 
fits into the local scale invariance scheme. 

In future work it will be interesting to extend our calculations to the response and 
correlation of Fourier modes of the defect density, which contain the global (energy) 
observable as the limiting case of zero wavevector. Such Fourier decompositions have 
proved valuable in the past in, e.g., Ising models [20] and also lend themselves naturally 
to a comparison with field-theoretical results in two or more dimensions [38]. Taking a 
broader view, it will be essential to establish how generic our finding of negative response 
functions is: the link to activated dynamics suggests that, once properly measured using, 
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e.g., global observables, negative responses may be much more common than previously 
thought. The corresponding negative FDRs, while not having a natural interpretation in 
terms of effective temperatures [5] , may be worthwhile in and of themselves in elucidating 
the nature of the activated dynamics. 



Appendix A. Special Functions 

Modified Bessel functions of integer order are most conveniently defined through the 
Fourier series of their generating function 

oo oo 

e tcos ^= /n(t)e i ^ = / (t) + 2^/ n (t)cos(^). (A.l) 

n=— oo n=l 

Note that useful resummation properties follow immediately. For instance, setting (p = 
gives e* = Io(t) + 2^^Li^n(£)- An explicit integral representation for the I n (t) is 
obtained by Fourier transformation, 

I n (t) = 1^ ^ cos(rup) e icos( ^. (A.2) 



2tt 

Elementary properties of the I n {t) apparent from these equations are I- n {t) = I n {t), 
In(~t) = (— l) n I n {t) and I n (0) = 5 n , t o. One shows that further 

d t I n (t) = ^[I n ^(t) + I n+1 (t)], (A.3) 

Th 1 

-I n (t) = -[I n ^(t)-I n+1 (t)}, (A.4) 

which we use extensively throughout the main text. From Eq. (lA.ip we obtain the 
convolution property 

oo 

I m (h + t 2 ) = Y In(tl)I m - n (t 2 ). (A.5) 



n=— oo 



The asymptotic behaviour of the modified Bessel functions of fixed and finite order n 
in the limit t — > oo is given by [57] 



e-*I B (f) ~ -L, (A.6) 
V znt 

whereas, if we simultaneously take t,n — > oo with n 2 /t fixed, 

e -%(t)~-^ e - n2 /^. (A.7) 

Another important function in our analysis is H n (t), a full characterization of which 
can be found in [51]; here we only summarize those features that are pertinent to the 
present context. First, H n (t) may be defined through 



1 



H n (t) = - I dre- T [I n -i(r) - I n+1 (r)} . (Ai 
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Clearly H_ n {f) = —H n (t) is odd in n and hence H (t) = 0. Remarkably, all functions 
H n (t) can be decomposed into modified Bessel functions via the recursion 

H n+1 (t) = H n (t) + 5 nfi - e~* [/„(*) + / n+ i(t)] , (A.9) 

which applies for n > and is derived by expressing H n+1 (t) — H n {t) via Eq. (IA.8I) 
and noting that the integrand becomes \ e~ T [— /„_i(r) + J n (r) + 7 n+1 (r) — / n +2( r )] = 
-d T e- T [I n (r) + I n+1 (t)]. Iterating Eq. ([Qjl we then have #i(t) = l-e - *[/o(t) + Ii(*)], 



i?2(t) = 1 — e~*[/o(i) + 2Ji(i) + i2(t)], etc. As a consequence i? n (t — ► oo) = 1 for all 
n > 1. The latter is useful for deriving asymptotic expansions of H n (t): we rewrite 
Eq. (IA.8I) as iJ„(t) = 1 — | J 00 dr e~ T [J n _i(r) — J n+ i(r)] so that r > t in the integrand. 
Then, using Eqs. (1A.41) . ( 1A.6I) and ( 1A.7I) we obtain that at fixed n > 1 and for t — > oo 

H n {t) ~ 1, (A. 10) 

while for t, n — > oo with n 2 /t fixed, 



# n (t) ~ $ ( -|=J where = -J= / dwe~ u \ (A.ll) 



This defines the symbol $(2) = 1 — erf(z) for the complementary error function. 
It remains to add that beyond Eqs. (1A.8I) and (IA.9I) there is yet another useful 
representation for the function H n (t), viz. 

H t2 ^(2t) = l- £}Gg(t), (A. 12) 

ii <j2 

which applies for zi < «2- Consistently with the notation elsewhere the sum in Eq. ( 1A.12I) 
runs over both G Z, but subject to the constraint j\ < j 2 - In order to establish 

the validity of Eq. ( 1A.12I) we substitute the Green's function Eq. (liSjl . set 32 = ji + m, 
shift j = ji — i\ and abbreviate n = i 2 — %\ > 1, which yields 



00 



if n (2t) = 1 - e- 2t [Ij(t)I j+m -n(t) ~ Ij-n^Ij+mit)]. 

m=l j=—oo 

After convolving the modified Bessel functions using Eq. (1A.5|) this expression reduces to 



H n (2t) = 1 — e 21 J2m=i [Im-n{2t) ~~ Im+n(2t)]- The remaining summation over m > 1 
is telescopic and one easily verifies that it obeys the recursion Eq. flA.9j) . thus proving 
the validity of Eq. flA~T2l . 

Appendix B. Subdominant Responses 

Here we establish subdominance of the corrections A^ s \t w ,t w ) to the instantaneous 
response functions for asymmetric and symmetric perturbations. These are given by 

+\ i+hh T S j>i2 -i) (elEi^m^m, e n \l). (B.l) 

The — and + alternatives for the signs correspond respectively to the asymmetric 
case, Eq. Q73]), and the symmetric one, Eq. (I8"TI) . To evaluate Eq. (IB.ip we substitute 
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hi — 1 — Ei t i + i as usual, which reduces the averages to combinations of the form 
Eqs. ( H0]HT1) . Then, making use of the recursion Eq. flA.91) yields 

Ag(tw,tw) = -~A id p i2 „ h (2t w ) and Ag(t w , t w ) = -B id Vi2 _ h (2t w ), (B.2) 

with and -Bij the coefficients defined below Eq. f!73|) and flHUl) . respectively, and 

<^ n (2t w ) = e~ 2t ™[I - I 2 ](2t w )H n (2t w ) + e- 4 *"[(Jo + - / n+ i)](2t w ). (B.3) 

At large t w and fixed n 2 /t w this scales like 

1 



<p n (2t^ 



2y/nf- 



3/2 



// \ A ii I n 2 



$ — = H = — = exp . 

1 2Vt^J ^2^/tZ F \ At 



(B.4) 



which reduces correctly to the scaling Lp n (2t w ) ~ l/(2y / 7rt w 3//2 ) for large t w but fixed 
n. Comparison to the leading contribution in the asymmetric response obtained from 
Eq. (!75|) shows that (t w , t w ) is uniformly - with respect to n = i 2 — i\ - subdominant 
and may therefore be ignored. The situation is somewhat different for the symmetric 
case: at fixed n and large t w the leading term in the symmetric response Eq. fl82|) scales 
like 0(n/t w 3 ^ 2 ) whereas A^ s j(t w ,t w ) = 0(l/t w 3 ^ 2 ). For small n > 1 contributions from 
A^j(t w , t w ) cannot be neglected, however, the leading term in Rfj(t w , t w ) grows linearly 
with n and thus quickly dominates A^ s j(t w , t w ). Specifically, once n 2 /t w = 0(1) we 
have Ag(t w ,t w ) = 0(l/tj' 2 ) while the leading term in Eq. ([E2D scales like 0(l/t w ). 
Altogether contributions from A^ s j(t w ,t w ) to the long time scaling of the symmetric 
instantaneous response become negligible when n — — i\ ^> 1. 
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